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The progress in our understanding of several aspects of turbulent Rayleigh-Benard convection is 
reviewed. The focus is on the question of how the Nusselt number and the Reynolds number 
depend on the Rayleigh number Ra and the Prandtl number Pr, and on how the thicknesses of 
the thermal and the kinetic boundary layers scale with Ra and Pr. Non-Oborbeck-Boussinosq 
effects and the dynamics of the large-scale convection-roll arc addressed as well. The review ends 
with a list of challenges for future research on the turbulent Rayleigh-Benard system. 



I. INTRODUCTION 

Rayleigh-Benard (RB) convection - the buoyancy 
driven flow of a fluid heated from below and cooled 
from above - is a classical problem in fluid dynam- 
ics. It played a crucial role in t he development of sta - 
bility theory in hydrod ynamics ( Chandrasekhad ( 198lh . 
iDrazin and R"ddl (|l98lh ) and had been paradigmatic in 
patter n formation and in the study of spatial-t emporal 
chaos (|Bodenschatz et al\ (|2000t ). lGetlh3 (|l998^ ). From 
an applied viewpoint, thermally driven flows are of ut- 
most importance. Exa mples are therma l conv ection in 
the atmosphere (see e.g . Hartmann et al. ' ('20011)), in the 
oceans (see e.g. [Marshall and Schott (1999)) (inclu ding 
thermohaline convection, see e.g. iRahmstorj (|200G[ )). m 
buildings (see e.g. iHunt and Lindeni " 1999( )). in process 
technology, or in metal-production processes (see e.g. 
iBrent ef all ( 1988t )). In the geophysical and astrophys- 
ical cont ext, we mention conve ction in Earth's mantle 
(see e.g. iMcKenzie aZI (Il974 )) in Earth's outer core 
(see e.g. ICardin and OlsonI (Il994l)). and in stars includ- 
ing our sun (see e.g. lCattaneo et all (200l)). Convection 
has been associated with the g eneration and reversal of 
Earth 's magnetic field, see e.g. iGlatzmaier and Robertd 
(1 19951) . 

Even if one restricts oneself to thermally driven flows 
in a closed box, there are so many aspects that not all 
of them can be addressed in this single review. We shall 
focus on developed turbulence when spatial coherence 
throughout the cell is lost and only on the large scale 
dynamics of the flow and aspects intimately connected 
with it such as the boundary layer structures. The scal- 
ing of the spectra of velocity and temperature fluctua- 
tions, or of the corresponding structure functions, will 
not be addre ssed . These issues had been discussed in 
the review bv lSiggial(ll994D . but meanwhile considerable 
progress has been achieved, in particular on the question 
of whether and w here in the flow to expect Bolgiano- 
Obuk ho y scaling ( Bolgianol ( 19591 ). iMonin and YaglomI 
( 19751 ) , Obukhov _j 1959!)) of t h e structure fu n ctions , 



ISun et al\ (|2006D . 

The question to be asked about the Rayleigh-Benard 
problem is as follows: For a given fluid in a closed con- 
tainer of height L heated from below and cooled from 
above, what are the flow properties inside the container 
and in addition, what is the heat transfer from bottom 
to top? Here spatially and temporally constant temper- 
atures are assumed at the bottom and top. In section 
mil we will discuss to what degre e this a ss umption can be 
i ustifi e d in reality fero wn et al\ ( 2005a| ). [Chaumat et al\ 
(I2OO2D , FVCTzicco (2004) ) . 

The problem is further simplified by t he so-called 
Oberb e ck-Boussinesq (OB) appro x imation (iBoussinesol 
(ll903D . lLandau and Lifshitd (ll987D . l0berbeckl (|1879D ) in 
which the fluid density p is assumed to depend linearly 
on the temperature. 



p(r) = p(To)(i~/3(T-ro)), 



(1) 



with [3 being the thermal expansion coefficient. In ad- 
dition, it is assumed that the material properties of the 
fluid such as /?, the viscosity i/, and the thermal diffu- 
sivity K do not depend on temperature. The govern- 
ing equations of the R B problem are then the Obe rbeck- 
Boussinesq equations ( Landau and Lifshitj ( 1987! )) 



dtUi -\- UjdjU 
dtO 



-dtP 



ujdj0 



(2) 
(3) 



J k rrr 



for the velocity field u{x,t), the kinematic pressure field 
p{x, t), and the temperature field 6{x, t) relative to some 
reference temperature. Here and in the following we as- 
sume summation over double indices; 5ij is the Kronecker 
symbol. The Oberbeck-Boussinesq equations are assisted 
by continuity diUi = and the boundary conditions 
It = for the velocities at all walls, 0{z = —L/2) = A/2 
for the temperature at the bottom plate, and 6{z = 
L/2) = —A/2 for the temperature at the top plate. At 
the side- walls the condition of no lateral heat flow is im- 
posed. The limitations of the Oberbeck-Boussinesq ap- 
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Within the OB approximation and for a given ccU ge- 
ometry, the system is determined by only two dimension- 
less control parameters, namely the Rayleigli number and 
the Prandtl number 



Ra 



Pr=^. 

K 



(4) 



The cell geometry is described by its symmetry and one 
or more aspect ratios F. For a cylindrical cell T = d/L, 
where d is the cell diameter. 

The key response of the system to the imposed Ra is 
the heat flux H from bottom to top. The dimension- 
less heat flux Nu = H/{AAL^^ is the Nusselt number. 
Here A = Cppn is the thermal conductivity. Within the 
Oberbeck-Boussinesq approximation one obtains for in- 
compressible flow 



Nu : 



(5) 



Here (.)^ denotes the average over (any) horizontal plane 
and over time. Correspondingly, {.)y used below denotes 
the volume and time average. 

Another key system response is the extent of turbu- 
lence, best expressed in terms of a characteristic velocity 
amplitude nondimensionalized by i/L^^ to define a 
Reynolds number 



Re = 



U 



vL- 



(6) 



As we shall see in section IIVI there are various reason- 
able possibilities to choose a velocity, e.g. the components 
or the magnitude of the velocity field at different posi- 
tions, local or averaged amplitudes, turnover times or 
frequency peaks in the thermal spectrum, etc. In some 
parameter ranges these amplitudes differ not only in mag- 
nitude but even show different dependences on Ra and Pr 
(jBrown et aLl (|2007^ . [Sugivama et aZ.I (|2008[ )). Mostly we 
shall restrict ourselves to that Reynolds number which is 
associated with the larg e-scale fiow (LSC), als o called the 
"wind of turbulence" U (iNiemela etU] (|200lh , ISun et al\ 
(l2005c[) . IXia et al\ (|2003f )). There is discussion in the lit- 
erature whether the LSC evolves out of the well-known 
cellular structure s at small Ra or whether i t does not. 
On the one hand iKrishnamurti and Howard! (Il98l[ ) per- 
formed experiments from which they concluded that the 
LSC is not a simple reminder and continuation of the 
roll structure observed just after the onset of convection. 
On the other hand we are not aware that their obser- 
vations have b een confirmed. Even explicit sea rch for 
such a mode as IKrishnamurti and" Howard! (|l98lh report 
was n ot su ccessful, see the sum marizing review bv lBuss^ 
(|2003[ ) and lHartlep et al\ (|2005[ ). These authors conclude 
that the LSC at large Ra indeed is a reminder of the low 
Ra structures. The dynamics of the large scale wind, its 
azimuthal oscillation, diffusion, reorientation, cessation, 
and possible breakdown at very large Ra will be discussed 



The key question to ask is: How do Nu and Re 
depend on Ra and Pr? The experimental situation 
will be the subject of sections IIHI for Nu(Ra,Pr) and 
rvl for Re(Ra,Pr). Results from numerical simula- 
tions are reported in section |Vl However, first (sec- 
tion ini we will summarize older theories (|H.A[) and 
then (subsection IH.BD the Grossmann-Lohse (GL) the- 
ory ([CrossmamandLohsi ((IMO, [1^ In 
subsection III. CI we will discuss theories about a possi- 
ble asymptotic regime at very large Ra and strict upper 
bounds for Nu. 

Section IVII is devoted to the structure and width of 
the thermal and kinetic boundary layers (BLs). The 
thermal BLs play a crucial role in determining the heat 
transfer, and the kinetic BLs provide viscous dissipa- 
tion of the L SC. Another important feature is the ther- 
mal p l umes (|Funfschilling and Ahlerd 12004). Kadanofl 
([200iyXi aZ.I(l2004D . IZhou et aLl (l2007bf ). IZocchi et al. 
(|l990l )) that detach from the thermal boundary layers; 
they contribute to the driving of the flow. In order to 
give an idea of the importance and organization of these 
plumes and their shapes (for large Pr) we show their 
shadowgra ph visualizat i on in the top image of figure [1] 
taken from lShang et al] (|2003[ ). The middle image of fig- 
ure [T] shows a streak picture of the temperature distribu- 
tion close to the upper plate, including a detaching plume 
for medium Pr. The lower image is a velocity-vector map 
of the LSC. 

As already mentioned, section IVIII is devoted to non- 
Oberbeck-Boussinesq effects and section IVIIII deals with 
the global wind dynamics. In section HXl we outline some 
major issues in Rayleigh-Benard convection for future re- 
search. 



II. THEORIES OF GLOBAL PROPERTIES: THE 
NUSSELT AND REYNOLDS NUMBER 

A. Older theories for Nu(Ra,Pr) and Re(Ra,Pr) 

For a detailed di scussion of th e theories developed prior 
to the review by Siggia (Il994l) we refer to that paper 
and to IChandrasekhaii (|l98lD 1 
power laws 



These theories predicted 



Nu - Ra'^''"Pr°'" 
Re - Ra'^^'Pr""^- 



(7) 
(8) 



for the dependences of Nu and Re on Ra and Pr. A sum- 
mary of predicted exponents is given in table IH Early ex- 
periments were of limited precision, and were consistent 
with power-law dependences over their limited ranges of 
Ra and Pr. 

The conceptually easiest early theory is Malkus' 
marginal-stability theory of 1954. It assumed that the 
thermal BL thickness adjusts itself so as to yield a criti- 
cal BL Rayleigh number. Th is immediately gave 'yNn = 
1/3. After the e xperiments bv lChu and GoldsteinI ( 19731 ). 
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Reference 


Pr & Ra range 




aNu 




"/Reflect 


CtRe 


Davis a922a,b~) 


Ra small 


1/4 










Malkus (1954) 




1/3 










Kraichnan (19621 


Ra ultimate, 
Pr< 0.15 


1/2 


1/2 


1/2 




-1/2 




Ra ultimate, 


1/2 


-1/4 


1/2 




-3/4 




0.15 < Pr < 1 












SDiegel (1971) 


Ra ultimate 


1/2 


1/2 


1/2 




-1/2 


CastainK et al. (1989) 




2/7 




1/2 


3/7 




Shraiman and Sieecia (1990) 


Pr> 1 


2/7 


-1/7 


3/7 




-5/7 


Yakhot (1992) 




5/19 






8/19 




Zaleski (1998) 




2/7 










Cioni et al. (1997) 


Pr< 1 


2/7 


2/7 




3/7 


-4/7 



T ABLE I Pow er-law exponents for Nu and Re as functions of Ra and Pr predicted by theories developed prior to the review 
bv lSiggial (|l994l ). The exponents are defined by Eqs. (O and ([8}. Whereas Re is based on the large-scale wind- velocity, Rcfiuct 
is based on the velocity fluctuations. 



experiments i n cryogenic helium (ICastaing al. (Il989f). 
iHeslot et all (Il987l). iProcaccia et al\ (|l99l[ ). ISano et al\ 
(ll989D . IWu et oi.l(ll99nh ) had suggested a smaller power- 
law exponen t, the Chic ago group developed the mixing- 
zone model dCastaing a/. (198lt) ) which later was ex- 
tended by ICioni et al\ ( 199'n ) to include the Prandtl- 
number dependences. The central result was jnu = 2/7. 
The same scaling exponent could also be obt ained from 
the BL theory of IShraiman and Siggial ( 1990l ). assuming 
a turbulent boundary layer. The assumptions of that the- 
ory are however very different from those of the mixing- 
layer theory, leading to very different power-law expo- 
nents for the dependences on the Prandtl number, see 
table H 

As we will see later, the assumption of a fully developed 
turbulent BL is far from being fulfilled in the parameter 
regime of the Chicago experiments. That can already be 
seen from an estimate of th e coherence leng t h £ of the RB 
flow. Taking the data from lProcaccia et al. (|l99l[ ) for the 
scaling of the velocity fluctuatio ns and of the crossover 
freque ncy to the viscous subrange. iGrossmann and Lohsj 
(I1993D obtained £/L ^ 50Ra~°-^^ For the F = 1/2 cell of 
IProcaccia et d \ (fl99lh this implies that only at Ra w 10 
the coherence length becomes about 1/3 of the lateral cell 
width and 1/6 of its height, a pre-requisite for indepen- 
dent fluctuations to develop in the bulk. Estimates based 
on £ « IOt], where rj is the (locally or globally defined) 
Kolmogorov scale, give similar results. The transition to 
turbulence in the BL is correspondingly expected only 
at much large Ra, namely, at Ra w 10^'^ (at the edge of 
the achievable regime in the Chicago experiments), as we 
will see in the next subsection. 

In any case, at "large enough" Rayleigh number a tran- 
sition should occur towards an ultimate Raylei gh- number 
regim e. Such a regime was first suggested by Kraichnaiil 
(| 19621 ). ISpiegell (|l97l[ ) hypothesized that in that regime 
the heat flux and the turbulence intensity are indepen- 



sivity, which leads to = 1/2 (for more details, see 
section |ILC[) . Though in those days (1971 and before) 
no measured power-law exponent was even close to that 
value, that paper has been extremely influential, perhaps 
also because from a mathematical point of view no lower 
strict upper bound th an = 1/2 could be proven to 
exist for finite Pr (see iDoering and Constanti J (|l996l )). 

As will be shown in sections IIIII and llVi the experi- 
ments of the last decade reveal the limitations of most of 
these older theories. 



B. Grossmann-Lohse theory for Nu(Ra,Pr) and Re(Ra,Pr) 

Given the increasing richness and precision of exper- 
imental and numerical data for Nu(Ra,Pr) (section WU^ 
and Re(Ra,Pr) (section HV| . it became clear near the end 
of the last decade that none of the theories developed 
up to then could offer a unifying view, accounting for all 
data. In partic ular, the pre d icted P randtl-number depen- 
dence s of Nu ( Cioni etall ( 1997t ). IShraiman and Siggial 
(IT990l) ) are in disagreement with measured and cal- 
culated data. The refore, i n a s e ries of four papers, 
IGrossmann and LohsG (,2000| , l200ll . l2002l I2004D tried to 
develop a unifying theory to account for Nu(Ra,Pr) and 
Re(Ra,Pr) over wide parameter ranges. 

The backbone of the theory is a set of two exact rela- 
tions for the kinetic and thermal energy-dissipation rates 
e„ and eg respectively, namely 



A2 



e„ = (i^(a,?/,(a^,t))2)^^ = — (Nu- l)RaPr-2, (9) 

(10) 



ee = {K{di6{x,t)f)^ ^ k—'Nvl 



These relations can easily be derived from the Boussinesq 
equatio ns and the corresponding boundary conditions 
(see e.g. lShraiman and Siggial ( 1990[ )). assuming only sta- 
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FIG. 1 Upper: Shadowgraph visualization of rising and 
falhng plumes at Ra = 6 .8 • 10^, Pr = 596 (dipropylene gly- 
col) in a r = 1 cell. From lShang et al\ lj2003r ). Middle: Streak 
picture of temperature sensitive liquid crystal spheres taken 
near the top cold surface in a F = 1 cell at Ra = 2.6 • 10^ 
and Pr — 5.4 (water), in order to visualize plum e detachment. 
The v iew shows an area of 6.5cm by 4cm. From lOu and TongI 
Bottom: Time-averaged velocity-vector map in the 

lane o f the LSC at Ra = 7.0 x 10^ (adapted from lSun et al\ 

20053)). 



is to split the volume averages of both the kinetic and the 
thermal dissipation rate into respective bulk and bound- 
ary layer (or rather boundary layer-like) contributions, 

e« = ^u.BL + £u,bulk: (11) 

ee = ee.BL + ^eMuik- (12) 
The motivation for this splitting is that the physics of the 




FIG. 2 Sketch of the splitting of the kinetic (upper) and 
thermal (lower) dissipation rates on which the GL theory is 
based. In both figures the large-scale convection-roll with typ- 
ical velocity amplitude U is sketched. The typical width of the 
kinetic BL is A^, whereas the typical thermal BL thicknesses 
and the plume thicknesses are Xg. Outside the BL/plume 
region is the background fiow "bg". 



pation rates is fundamentally different and thus the cor- 
responding dissipation rate contributions must be mod- 
eled in different ways. The phrase "BL-like" indicates 
that from a scaling point of view we consider the detach- 
ing thermal plumes as parts of the thermal BLs. Thus 
instead of BL and bulk we could also use the labels pi 
(plume) and bg (background) for the two parts of the 
dissipation rates. A sketch of the splitting is shown in 
figure [21 Rather than eq. one therefore could also 
write 



(13) 



signalling the contributions from the BL and the plumes 
(pi) on the one hand and from the background (bg) on 
the other. 

Two further assumptions of the GL theory are indi- 
cated as well in figure [U namely that there exists a large- 
scale wind with only one typical velocity scale U (defin- 
ing a Reynolds number Re = UL/v), and that the kinetic 
BLs are (scalingwise) characterized by a single effective 
thickness A„ regardless of the position along the plates 



the velocity scales and in section I VII for the BL thick- 
nesses, both assumptions are simplifications. In particu- 
lar, even the scaling of the kinetic BL thickness with Ra 
may be different at the side-walls as compa r ed to the top 
and bott om plates (se e lLui and Xial (Il998l). ^u and Xial 
(ll998bD . lxi^ and Xial (1X9971 ). IXin et al\ (|1996D ). Never- 
theless, for the sake of simplicity and in view of Occam's 
razor, - and consistent with the rnost rec ent experimental 
results for the BLs bv lSun et all ( 2008f ) - these simplifi- 
cations have been used. 

Accepting the splitting ([TT]) and HI]) (or ((T3])), the im- 
mediate consequence is that there are four main regimes 
in parameter space: Regime I in which both e„ and eg are 
dominated by the BL/plume contribution, regime II in 
which e„ is dominated by eu,buik and ee by ee,BL, regime 
III in which is dominated by eu.BL and ee by eeMk, 
and finally regime IV in which both and eg are dom- 
inated by their bulk contributions. It remains to be de- 
termined where in Ra-Pr parameter space the crossovers 
between the different regimes are located. 

The next step is to model the individual contributions 
to the dissipation rates. We start with the bulk contribu- 
tions. The turbulence in the bulk is driven by the large- 
scale wind U. The corresponding time scale therefore 
is L/U, and from Ko lmogorov's energy-cascade picture 
(see e.g. iFrischI (|l995[ )) the bulk energy dissipation-rate 
scalingwise becomes^ 



^u.bulk 



El 

L 



(14) 



This seems justified because the turbulence i n the bulk 
is mo re or less homogeneo us and isotropic ( Sun et oil 
(|2006l ). IZhou et al\ (|2008al) V The same reasoni ng can 



be app lied to the temperature equation, see again iFrischI 
( 19951 ). The bulk thermal dissipation-rate then becomes 



A2 

— -PrRe. 



(15) 



The scaling of the boundary-layer contributions to the 
dissipation rates are estimated from their definitions 
as BL averages e„,_BL ~ {{diUj{x £ BL,t))^')y and 
£0,BL = K{{di9{x e BL,t)Y)y, namely 



Am 



and 



L ■ 



(16) 



(17) 



As detailed by iGrossmann and Lohsd (I2004D . the ki- 
netic and thermal BL thicknesses A„ and Xg are ob- 
tained from the Prandtl-Blasius BL theory ( BlasiusI 



(I1908D. ICowlevI (l200lh. i Meksvnl (Il96l[ ). iPrandtl (I1905D . 
ISchhchting and GerstenI (|200(3l )): 



Au 

T 



(18) 



where a is a dimensionless prefactor of order one, and 
A« 



L 



i?e~^/^Pr~^/2 for Pr < 1, (19) 
i?e"^/2p^-i/3 forPr>l. (20) 



Note that scalingwise laminar BL theory is applied which 
seems justified because of the low prevailing boundary 
Reynolds numbers. Further below it will be estimated 
when this assumption breaks down for increasing Re. 
In the small Pr regimes (eq. (|19p ) (label "1" stands for 
"lower" in Fig. [3]) the kinetic BL is nested in the thermal 
one, A„ < Ae, whereas in the large Pr regimes (eq. (PO]) ) 
("u" for "upper" in Fig. [3]) the thermal BL is nested 
in the kinetic one, \g < A„. The transition from one 
regime to the other is modelled "by hand" through a 
crossover function f{xg) = {1 + xf)~^/^ of the variable 
xg_=\u/\g = 2aNu/Re^^^, see IGrossmann and Lohsd 
( 200ll ). Note that in the crossover function f{Xu/Xg) the 
thermal BL thickness Xg has been replaced by L/(2Nu). 
Finally, when Re becomes very small the expression (jl8p 
for the kinetic BL thickness diverges, while the physical 
A„ is limited instead by an outer length scale of the or- 
der of the cell height L. This saturation is happening at 
some small but a priori unknown Reynolds number Rec- 
The transition towards the saturation regime is again 
modeled "by hand" with the crossover function g{xL) = 
XLil + xi )-^/^ with XL = XJRe )/XJRer) = ^RejRe 



(again see Grossmann and Lohsd ( 200l[ ) for details). 

When putting the splitting and modeling assumptions 
together with the two exact relations ^ and (fTO|) . one 
finally obtains two implicit equations for Nu(Ra,Pr) and 
Re(Ra,Pr) with six free parameters a, i?ec, and Ci, i ~ 
1,2,3,4: 



{Nu - l)RaPr-^ = ci 



Re^ 



gi^RejRe) 



^ csRe^/^Pr^'^ 

^ ( 2aNu 
+ CiPrRef —j=g 



\\/Rec \ V Re 




C2i?e^ (21) 



-,1/2 



(22) 



The —1 on the Ihs of eq. (|22p stems from the contri- 
bution of the molecular transport, which survives when 
the Peclet number Pe = RePr = UL/k tends to zero, 
Pe 0, cf. [G rossmann and Lohsd (|2008l ). This happens 
if either the velocity field decreases, Ui — 0, or if the ther- 
mal diffusivity becomes large, k oo. In either case the 
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form dj9 = 0, whose solution with the proper boundary 
conditions is 9 = —AL^^z. Inserting this sohition into 
the Nusselt number definition ([5]) gives hmpe_,o Nu ~ 1. 
Of course the —1 does not matter much in the turbulent 
regime with large Nu. 

The six parameters in eqs. (|2ip and (|22p were adjusted 
so a s to provide a f i t to 1 55 data points for Nu(Ra,Pr) 
from lAhlers and These data were in the range 

3-10^ < Ra < S-IO^ and 4 < Pr < 34 for a T = 1 cylindr i- 
cal cell. As elaborated bv lGrossmann and Lohsd ( 2002h . 
in order to fix the parameter a o ne also needs to know 
Re for (at least) one pair (R a,Pr); Grossmaiin and Lohsd 
(|2002l ) took that value from lQiu and Tond l|200lU ). The 
final results were a = 0.482, ci = 8.7, C2 = 1.45, 



C3 



0.46, 



C4 



0.013, and Rcr = 1.0. With this set 



the data of lAhlers and Xul (|200l[ ) were described very 
well. Later these data were adjusted for side-wall and 
plate corrections. However, t he agreement with t hem a s 
well as with additional data ( Funfsehilling et oil (j2005[ )') 
for Ra up to 3 • 10^" and Pr = 4 .38 (see figure S]) is 
still excellent (see subsection IIH . C P . For the Nu(Ra,Pr) 
and Re(Ra,Pr) predicted with these parameters over wide 
ranges of Ra and Pr we re fer t o the figures given by 
iGrossmann and Lohsd (l200lh and[G rossmann and Lohsd 
( 20021 ). We mention that in principle one expects an 
aspect-ratio dependence of the Ci, since the relative con- 
tributions of BL and bulk change with aspect ratio. How- 
ever, from experiment it is known that the F dependence 
of Nu(Ra, Pr) is very weak in the explored range of Ra 
and Pr (see subsection IHLEp . 

After the determination of these six parameters, 
Nu(Ra,Pr) and Re(Ra,Pr) are given for all Ra and Pr 
by equations (PT|) and (I22p . In addition, also the Ra- 
Pr parameter-space structure with all transitions from 
one regime to another is determined. The corresponding 
phase diagram is reproduced in figure [3] 

One central assumption of the GL theory is the ap- 
plicability of the scaling of the Prandtl-Blasius lami- 
nar BL theory. For increasing Ra and thus increas- 
ing Re this assumption will ultimately break down; 
the BLs are expected to becom e turbulent as well. 
IGrossmann and Lohsd (|2000l . I2002D provide an estimate 
for the Rayleigh number at which the breakdown occurs, 
based on the shear Reynolds number Rcg = XuU/v ~ 
lie. For RB experiments using classical fluids over 



"typical" Ra and Pr ranges Re^ is not particularly large. 
This reflects the relatively low degree of turbulence in the 
interior, which also beco mes evident from flow visualiza- 
tions li ke those d o ne by iTilgner et al. (Il993l). Xi et al. 



20041) . IXia et al\ (|2003h . and iFunfschiUing and Ahleri 
20041 ). For example, with Pr ~ 4 one has Rcg = 15 when 
Ra = 10^ and Re w 900, and Re,, = 190 for Ra ^ 10^^ 
and Re ~ 140000. The dashed line in figure [3] is based on 
the critical value i?e* ~ 420. Beyond i?e* the kinetic BLs 
become fully turbulent and the Prandtl-Blasius scaling is 
no longer applicable. It is not totally clear what will hap- 
pen in this ultimate regime of thermal convection. That 



Note that the notion of laminar kinetic BLs in RB flow 
should not be confused with time independence or lack of 
chaotic behavior. The detaching thermal plumes intro- 
duce time dependences an d chaotic behavior into the k i- 
netic BL; but, as shown bv lGrossmann and Lohsd ( 20041 ). 
also then the Prandtl-Blasius scaling laws for the thick- 
nesses of the BLs still hold. Assu ming a turbulent BL 
alread y for Ra ^ 10^^ as done bv IShraiman and Siggial 
(|1990D ' eads to a dependence of Nu on Pr that disagrees 
with experiments and numerical simulations. 

A detailed comparison of the GL theory with various 
data will be done in sections IIIII and IIVI Here we only 
stress that the theory has predictive power: The deter- 
mination of the free parameters was done in the limited 
parameter range 3 • 10^ < Ra < 3 • 10^ and 4 < Pr < 34, 
see blue stars in figure [31 The predictions of the theory, 
however, hold over a much larger domain in the Ra-Pr 
parameter space. 

We further note that thanks to eq. ^ the knowledge of 
the Nusselt number allows for an estimate of the volume 
averaged energy dissipation rate and derived quantities. 
E.g., when taking th e conditions of t he Oreg on cryogenic 
helium experiment ( Niemela et al\ ( 2000al )). for Ra = 
10^° and Pr = 0.7 one obtains either directly from exper- 
iment or from the GL theory a Nusselt number of 120 and 
with the experimental values L=lm and i/ ~ 5-lO^^m^/s 
an energy dissipation rate of = 3 • 10^*m^ / . At 
Ra = lO^'' and Pr = 0.7 one obtains Nu w 2400 and 
with = 10~^m^/s a value of e„ = 5 • 10~^m^/s'^. Both 
of these energy dissipation rates are about three orders 
of magnitude smaller than in typical wind tunnel experi- 
ments. From the volume averaged energy dissipation rate 
eq. ^ one can also obtain global estimates for the spa- 
tial coherence length £ which typically is about 10 times 



the Kolmogorov length scale i] 



.3/4/1/4 



/ej . E.g., for 



cryogenic helium (Pr = 0.7) at Ra = 10^ one obtains 

e/L w lOrj/L = lOPr^^VK^" ^ ly/'^Ra^^'^] w 0.08, 
which is small enough to allow for the loss of spatial 
coherence and the onset of turbulence in the bulk. In 
contrast, for the same Ra in water (at Pr = 4) one 
has e/L « 0.18 and in glycerol (at Pr=2000) even at 
Ra = 10^ one only has t/L « 0.9, so that there is no 
developed turbulence. In glycerol, only at Ra = 10^^ 
one obtains ijh < 0.2 and thus developed turbulence, 
according to this GL-model based estimate. 

Finally, we note that the GL approach also has 
been applied t o oth e r geom etries and flows: E.g. 
lEckhardt etd\ (|2000l l2007allbl) app l ied it to Taylor - 
Couette and pipe flow and lTsai oZ] (|2007l l2003l . l2005l ) 
to turbulent electroconvection. 



C. Is there an asymptotic regime for large Ra? - and strict 
upper bounds 



iKraichnanI (|l962l ) and later ISpiegell (|l97ll ) postulated 
an "ultimate" , or asymptotic, regime in which heat trans- 
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regime 


dominance of 


BLs 


Nu 


Re 


/; 

lu 

/< 
/> 

OO 


^u,BL, (-e.BL 


Au < Xe 

Au > Xg 

Xu = L/4 > Xe 
Xu = L/4 > As 


Ra'^' 
Ra'^^ 


Ra^/'^Pr-'^/^ 

^„l/2p^-5/6 

i?a2/3p,.-i 
Ra^/^PT-^ 


III 

IIu 


iu.bulk, ^e,BL 


Xu < Xg 
Xu > Xg 


Ra^^'Pr^'" 
Ra^'^ 


Ra^/'^Pr-^/'^ 
Ra^^Pr-^/s 


IIIu 
IIIoo 


tu,BL, ^e.bulk 


Xu > Xg 

Xu = L/4 > As 


Ra^/'' Pr-^/'' 
Ra^/'' 


Ra'^/rpr-e/7 
Ro^/^Pt-^ 


IVi 

IVu 


^Uybulkj ^0 ,bulk 


Xu < Xg 
Xu > Xg 


Ra'^'^Pr^^^ 
Ra^/'' 


^„4/9p^-2/3 



TABLE II The pure power laws for Nu and Re in the various regimes. Table taken from lGrossmann and Lohsd (|200ll L 



the kinematic viscosity and the thermal diffusivity. The 
physics of this ultimate regime is that the thermal and 
kinetic boundary layers, and thus the kinematic viscosity 
v and the thermal diffusivity k, do not play an explicit 
role any more for the heat flux. The flow then is bulk 
dominated. With proper non-dimensionalization, and in- 
cluding logarithmic corrections due to viscous sub-layers 
induced by no-slip boundary conditions, Kraichnan's pre- 
dictions for this regime read 

iVit - i?a^/2 (log i?a) "3/2 p^i/2^ (23) 
i?e - i?ai/2(logi?a)-i/2 Pr-i/^ (24) 
for Pr < 0.15, while for 0.15 < Pr < 1 he suggested 

Nu - i?ai/2(log i?a) '3/2 p^- 1/4^ (25) 

(26) 



Re - i?ai/2 (log i?a)- 1/2 Pr" 



3/4 



The Ra-number dependences agree with the depen- 
dences in regimes IV ; and IV, of the GL theory 
(iGrossmann and Lohsd (|2000L l200lL l2002l I2004D L ex- 
cept for the logarithmic corrections. The Pr-dependence 
within the GL theory in the ultimate regimes IV; and IVJ 
is different: 



Nu - Ra^/^Pr^/^, 



Re - Ra^/'^Pr-^/'^. 



(27) 



(28) 



Equation (gT]) was derived first bv lSpiegell (|l97lD from a 
model for thermal convection in stars. 

To illustrate the physical implications of the existence 
of the ultimate regime, Andreas Acrivos (priv. commu- 
nication) suggested the following Gedankenexperiment: 
Consider RB convection in a very large aspect ratio sam- 
ple, with the lateral dimension (say, the diameter D of 
a cylinder) much larger than the sample height L. Now 
fix all dimensional parameters (A, k, z/, /?, g, and D) 
and increase the sample height L, starting from zero, but 
such that always still D ^ L, i.e., remain in the large 
aspect ratio limit. How does the dimensional heat flux 



to Nu — 1. With increasing L, the decrease will be- 
come weaker. The regime Nu ~ Ra^/^ corresponds to 
the dimensional heat flux H being independent of L. For 
even further increase of L (with still D ^ L), the exis- 
fence of the ultimate regime Nu ^ Ra ' would imply 
that the dimensional heat flux H would increase again, 
namely with L^/^. This feature may be considered as 
counter-intuitive, however, our interpretation of the ul- 
timate regime (if it exists) is that with fully developed 
turbulence in the bulk, the increasing sample height L 
allows for larger and larger eddies which thus can trans- 
port more and more heat from the bottom to the top 
plate. 

Ever since Kraichnan's prediction in 1962, researchers 
have tried to find evidence for this regime. Various ex- 
perimental efforts will be discussed in subsection IIII.GI 

There are also numerical indications of the ultimate 
regime: In order to obtain a Nu ^ Ra^^^ power law, the 
classical velocity and temperature bo undary-conditions 
of the RB problem have been modifi ed: iLohse and Toschl 
( 20031 ) and lCalzavarini et al. I (|2005D performed numerical 
simulations for so-called "homogeneous" RB turbulence, 
in which the top- and bottom-temperature boundary- 
conditions have been replaced by periodic ones, with an 
unstratified temperature gradient imposed. The idea was 
to eliminate the BLs in th is way. The numerical results 
of lCalzavarini et all (|2005[ ) - including the found Prandtl 
number dependence ~ are consistent with the ultimate 
scaling eqs. (P7)) and (pS)) . where the Reynolds number is 
th at of the velocity fl u ctuati ons. As has been pointed out 
by ICalzavarini et all ( 20061 ). however, one should note 
that the dynamical equations of homogeneous RB turbu- 
lence allow for exponentially growing (in time) solutions, 
i.e., homogeneous RB turbulence does not have any strict 
upper bound for Nu. 

Such upper bounds do exist for the classical 
RB problem. Build ing on H oward 's sem inal vari- 
ational formulation (iHowardI (|l963l ll972D 'l. iBussd 
(|l969l ) coul d prove that Nu < (Ra / 1035V /^ for any 
Pr. Later iDoering and ConstantinI (|l996( ) derived a 
strict upper bound given by Nu < 0.167Ra^/^ — 
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( Doering and ConstantinI ( 1992| )). The hitherto abso- 
lute best asymptotic upper bound on Nu(Ra) comes 
from Plasting and KersweU ( 2 003.) , obtaining Nu < 
1 + 0.02634i?ai/^ which is 20% lower than Busse's 
best estimate. For arbritary Pr no power-law expo- 
nent of Ra smaller than 1/2 could hitherto be ob- 
tained as an upper bound. However, for infinite Pr 
IConstantin and Doerind ( 1999f ) could prove that Nu < 
const X Ra^/'^llogKa]^^^ . This res ult was improved late r 
to Nu < 0.64 4 xi^Q^/^MRal^/^ bv lDoering et aLl(|2006f) . 
lOtero et al\ (|2002t) obtained a strict upper bound for Nu 
for RB convection with constant heat flux through the 
plates (rather than with constant temperatures of the 
plates), namely, Nu < constxRa^^^ also for this case. We 
note that the scaling laws resulting from the GL theory 
are compatible with the upper bounds, including those 
in the large- Pr limit. 



system. Another unique property of materials at low 
temperatures is the extremely small heat capacity and 
large thermal diffusivity of the confining top and bot- 
tom plates which permit the study of temperature fluc- 
tuations at the fluid-solid interface when the heat cur- 
rent is held constant and led to the observation o f chaos 
i n a system governed by continuum equations ( AhlersI 
(|1974 |1975[ )V Additional advances m recent times 
have been due to the application of precision measure- 
men ts, using classical liquids a nd ga se s at ever increasing 
Ra (iFleischer and GoldsteinI (120021) . iFunfschilling et al\ 



(I2005D. iNikolaenko al\ (120051). iRoche et al\ (|2002f ). 
ISun et all ( l2005al). IXu e t al) ('200 01)) and oy er a wide 
range of Pr (|Ahlers a~Xu (,200lf )7 lxia et al\ (|2002f )). 



B. Side-wall and top- and bottom-plate-conductlvlty 
effects on Nu 



III. EXPERIMENTAL MEASUREMENTS OF THE 
NUSSELT NUMBER 

A. Overview 

During the last two or three decades measurements 
of Nu(Ra) as a function of such parameters as F, 
the extent of departures from the OB approxima- 
tion, the deliberate suppression of the large-scale cir- 
culation (LSC) by internal obstructions, the rough- 
ness of the confining solid surfaces, or deliberate mis- 
alignment relative to gravity have revealed various as- 
pects of the heat-transport mechanisms involved in 
this system. These efforts received a significant boost 
when it was appreciated that liquid or gaseous he- 
lium at low temperatures offered experimental oppor- 
tunitics not a vailable a t ambient temperatures ("Ahlers 
1974. igTSl . Behringer ( f 9851 , Nicmela and Srcenivasan 
2006b| ). frhrelfalll (|l975h ). Extensive low-temperature 
measurements of Nu(Ra) we re init iat ed by the 



Chica g o group (ICastaing et al\ ( 19891 ). iHeslot et al\ 
(Il987l). ISano et al\ (ll989^). followed by t he Grenoble 
group (IChavanne et al\ (Il997l Il996l . l200lh . iRoche et al] 



2001bl. I2OO2I I2004D') and t he Oregon/Trieste group 
(Niemcla et al. ' ('2000al . I2OOII). iN icmcla and SreenivasanI 
2003a . ,2006a , b). Niemela et ah (20001^,^). Among the 



advantages of the low-temperature environment is the ex- 
ceptionally small shear viscosity of helium gas which, at 
sufhciently high density, permits the attainment of ex- 
tremely large Ra. Further enhancements of the achiev- 
able Ra have be attained near the critical points of sev- 
eral fluids, including helium, where the thermal expan- 
sion coefficient diverges and the thermal diffusivity van- 
ishes, yielding a diverging Ra at constant A. Here, how- 
ever, it must be noted that on average the increase of 
Ra is accompanied by an increase of Pr (see Fig. [SI bot- 
tom) because Pr diverges as well at the critical point. 
This makes it difficult to disentangle any influence of 



A serious problem for quantitative m easurements o f 
Nu is the influence of the side wall (lAhlersI (l2000|) . 
Niemela and SreenivasanI ( 2003a[ ). iRoche et al\ ( 2001a ). 
Verziccd (|2002t) ). The wall is in thermal contact with the 
convecting fluid and shares with it, by virtue of the ther- 
mal BLs, a large vertical temperature gradient near the 
top and bottom and a much smaller gradient away from 
the plates. Thus, the current entering and leaving the 
wall is larger for the filled sample than it is for the empty 
one. In the wall near the top and bottom ends there is 
also a lateral gradient that will cause a part of the wall 
current to enter the fluid in the bottom half of the sample, 
and to leave it again in the to p half. This will influence 
the det ailed nature of the LSC ( Niemela and SreenivasanI 
(l2003aD ). however, the global heat current is determined 
primarily by processes within the top and bottom BLs. 
Thus it is insensitive to the detailed structure and inten- 
sity of the LSC and is not influenced much by this compli- 
cated lateral heat flow out of and into the wall. Therefore 
the problem reduces primarily to determining the current 
that actually enters the fluid. Approximate models that 
provide a corr e ction f or this wall e f fect ha ve been pro- 
posed (|Ahlersl (|2000D . iRoche et al\ ^Mk)^ but these 
are of limited reliability when the effect is large. The 
cryogenic measurements have a disadvantage because the 
sample usually is contained by steel side walls that have 
a relatively large conductivity A„, ~ 0.2 W/ (m K), while 
the fluid itself has an exceptionally small conductivity of 
order 0.01 W/(m K), giving A.u,/A ~ 20. In this case 
the models suggest that the correction is about 10% of 
Nu when Nu ~ 100 (Ra ~ 4 x 10^) and of course larger 
for smaller Nu. Even for Ra ~ 10^^ where Nu ~ 280 
a correction of about 6% is suggested. The net result 
is that the measured effective exponent of Nu(R a) is re- 
duced below its true value by about 0.02 or 0.03 (jAhlerd 
(I2OOOD ). ^'or gases near ambient temperatures with typ- 
ical thermal conductivities near 0.03 W/(m K), such as 
sulfur hexafluoride (SFg) and ethane (C2Hg), confined by 
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W/(m K) (|Ahi ers et al\ (|2007[ )). one approaches the case 
of perfectly conducting lateral boundaries where subtrac- 
tion of the current measured for the empty cell actually 
becomes a good approximation. Nonetheless results for 
Nu, although very precise, can not be expected to be 
very accurate. An exceptionally favorable case is that 
of wat er confined by relatively thin plastic walls ( Ahlerd 
( 2000l )). where Au,/A ~ 0.3. In that case the side-wall 
correction can be as small as a fraction of a percent and 
may safely be ignored for most purposes. An intermedi- 
ate case, for which reasonably reliable corrections can be 
made, is that of organic fluids confined b y various plastic 
walls w hich typically have A^/ A = 0(1) (jAhlers and Xul 
n the case of liquid metals, which are of interest 
because they have very small Prandtl numbers of order 
10~^ or less, A^/A is small (~ 2 for Hg and ~ 0.2 for Na 
as examples) and again the wall corrections are small or 
negligible. 

A second problem involves th e finite conduct i vity A p 



of the top and bottom plates (jChaumat et al\ (|2002t ). 
IVerziccol (|2004[) ). One would like X = ApL/(eANu) to be 
very large (here e is the thickness of one plate). Else the 
emission of a plume from the top (bottom) boundary will 
leave an excess (deficiency) of enthalpy in its former lo- 
cation, generating a relatively warm (cold) spot near the 
plate where the probability of the emission of the next 
plume is diminished until this thermal "hole" has diffused 
away by virtue of the plat e conductivity . This issue was 
explored experimentally bv lBrown et 'an (l2005aD bv mea- 
suring Nu(Ra) with high precision using water (A ~ 0.6 
W/(m K)) as the fluid and first Al and then Cu top and 
bottom plates of identical shape and size (see Fig. [5] be- 
low). The conductivities Ap^cu — 400 W/(m K) of Cu 
and Ap^Ai — 170 W/(m K) of Al differ by a factor of 
about 2.3 and thus yield different reductions of Nu(Ra) 
below the ideal value Nuoo for isothermal boundary con- 
ditions. The results permitted the extrapolation of Nu 
to Nuoo by the use of the empirical formula 

Nu = /(X)Nuoo; fix) = 1 - exp{-{aX)'') . (29) 

The parameters were found to be a = 0.275 and b = 0.39 
for L = 0.50 m, and f{X) was found to be closer to unity 
for smaller L. At fixed L both a and b (and thus f(X)) 
were independent of F. This plate-conductivity effect 
is expected to be relatively small for the cryogenic and 
room-temperature compressed-gas experiments because 
typically Ap/A = OiW^) and larger and thus X is very 
large unless Nu becomes extremely large. At modest Ra, 
say Ra < 3 x 10^, it is small also for Cu plates and organic 
fluids where Ap/A = 0(10'^). The plate correction is a 
very serious problem for measurements with liquid metals 
where for instance Ap/A ~ 50 for Hg and ~ 5 for Na. It 
has been suggested that this problem might be overcome 
by using a composite plate containing a volume partially 
filled with a liquid of high vapor pressure. In that case 
the condensation and vaporization of this fluid inside the 
plate can yield an effective plate conductivity much larger 



has not yet been implemented. 

The problem of the plate corrections was ad- 
dr essed again r ecent l y through numerical simulation s 
bv lAmati et al\ (|2005[ ) . l\^zicco and SreenivasanI (|2008f ). 
Results for Nu(Ra) obtained with constant heat-flux 
boundary conditions (BCs) at the lower plate and 
constant-temperature BCs at the upper plate were com- 
pared with Nu(Ra) for constant-temperature BCs at 
both plates. As shown below in Fig. [5^, the re- 
sults for both BCs agreed reasonably well with each 
other and with experiment up to Ra « 10^. This 
is also found by comparing 2D numerical simula- 
tions with c onstant temperature and constant flux BCs 
([Johnston and Doe ring (2007i ) and unpublished results 
by the same auth ors) . However, for Ra > 10^ 



(jAmati et al\ (|2005l ). IVerzicco and SreenivasanI (20081 )) 
the constant-temperature results for Nu were up to 20% 
higher than the constant heat-flux simulations and up 
to nearly 30% higher than experimental results. This 
suggests that in the high-Ra experiments, even after cor- 
rection for the finite plate-conductivity using Eq. (|29p . 
the results correspond more nearly to constant-flux at 
the bottom plate than to constant-temperature bound- 
ary conditions. An elucidation by experiment of this sur- 
prising result would be very welcome. 
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C. The Nusselt number for Pr ~ 4.38 obtained using water 
as the fluid 



For Pr ~ 4.4 and F = 1.00 l iiRh-a ccuracy measure- 
ments of Nu (|Funfschilling et al\ (|2005l )) for lO'^ < Ra < 
10^^ using water and copper plates are shown in Fig. 3] as 
circles. We focus on these data because for them the side- 
wall corrections are negligible and top- and bottom-plate 
corrections based on experiments with plates of different 
conductivities were made (see Sec. IIIFBp . The wide Ra 
range was achieved by using three samples with different 
L. The data before the plate correction are given as open 
circles. Corrected data are presented as solid circles. 

For F = 1 and Pr ~ 4 the experiment reaching the 
largest Ra was conducted using Cu plates and a wa- 
ter sample with L = 100 cm and reached Ra ~ 10^^ 
( Sun et These data are shown as open 

squares in the figure. It is gratifying that they are re- 
markably consistent with the open circles. However, the 
authors used an empirical plate correction with a = 0.987 
and b = 0.30 which yielded the solid squares in the fig- 
ure. In an attempt to develop an estimate of the un- 
certainty, we applied a correction using eq. ([^^ and the 
parameters a = 0.275 and b = 0.39 obtained from the 
L = 0.5 m sample. This yielded the up-pointing trian- 
gles. This correction is too small because measurements 
with a L = 0.25 m sampl e and the L = 0.50 m sam- 
ple bv lBrown et Ml (|2005aD revealed that the correction 
increases with L. Arbitrarily assuming a power-law de- 
pendence a = aoL^" and b = boL^'', an extrapolation to 
L = 1 m yielded a = 0.221 and b = 0.264, and via eq. ^ 
led to the down-pointing triangles. We consider the up- 
pointing and down-pointing triangles to be estimates of 
lower and upper bounds on the actual Nutxj. Arbitrar- 
ily adjusting a and b to the intermediate values 0.25 and 
0.32 respectively yielded the solid diamonds which are 
consistent with the data from the L = 0.5 m sample. 
New measurements in this very large cell with Al plates, 
which together with the Cu-plate data will yield better 
values of a and 6, are anxiously awaited. 

The solid line in Fig. [H is the GL prediction 
( Grossmann and Lohsd (|200l[ )). It gives the shape of the 
experimentally found Nu(Ra) very well for Ra < 10^". 
For larger Ra the data suggest 7e// = 1/3 whereas the 
model only reaches such a value for "feff as Ra oo 



where the model is no longer expected to be applicable. 
A detailed discussion of a number of other measure- 



ments for Pr = 0(1) and F = 


0.5 or 1 (lAhlers and Xu 


(2001 


). IChavanne et al\ (l200lh. Fleischer and Goldstein 


(2002 


).lNiemela et al. (2000a). 


Niemela and Sreenivasan 


(2003al). Nikolaenko and Ahlers 


(20031). Nikolaenko et al. 


(2005). iRoche et al\ (200l l2004^. Sun et al\ (2005a^. 


Xu et al. ( 20001)1 is beyond the scope of this re- 



view, although we shall re-visit a few of them be- 
low in Sec. IIII.GI We refe r to pr ev ious publications 



in the range Ra < 10^^ there are differences of up to 
20% or so between some of them. It is not clear whether 
the origin of these differences is to be found in exper- 
imental uncertainties, perhaps associated with wall or 
plate corrections o r other experimental effect s , or, as has 
been suggested by iNiemela and Sreenivasa 3 (l2003aD . m 
genuine differences of the fluid dynamics of the various 
samples. We find the latter explanation somewhat un- 
likely because, as discussed below in Sec. IIII.Fl the heat 
transport is determined primarily by boundary layer in- 
stabilities and is relatively insensitive to the structure of 
the LSC. 



D. The Prandtl-number dependence of the Nusselt number 

Fluids with Pr > 1 are plentiful in the form of various 
liquids, although accurate determinations of Nu(Ra) are 
in many cases problematic because the required physical 
properties are not known well enough. Typical gases not 
too close to the critical point have Pr = 0(1). The range 
Pr < 0.7 is difficult to access because most ordinary flu- 
ids have Pr greater than or c lose to the hard-sphere-g as 
value 2/3 (see, for instance, iHirschfelder et aU ( 19641 )'). 
Liquid metals, by virtue of the electronic contribution to 
the thermal conductivity, have Pr = 0(10"^) or smaller, 
leaving a wide gap in the range from 10~^ to 0.7. For 
the liquid metals it is difficult to obtain very large val- 
ues of Ra because the large thermal conductivity requires 
very large heat currents and tends to yield small Rayleigh 
numbers unless very large samples are constructed. An- 
other problem for liquid metals (see Sec. IIII.Bp is the 
uncertainty introduced by a large plate correction; how- 
ever, side-wall corrections should be negligible. 

In spite of these difficulties, several researchers at- 
tempted low-Pr measurements of Nu, in order to study 
the Pr dependence. Measurements w ith mercury (Pr = 
0.025) we re done bv iRossbvl (1 19691) (2 ■ 10^ < Ra < 



5-105), bv lTakeshita eraZ.I (Il996l) andlNaert al. 



ail dl995l 



199i 



(10^ < Ra < 10^), bv ICioni et 

(5 • 10^ < Ra < 2 • i qQ), and by I Glazier et al\ , . , 

(2 • 10^ < Ra < 8 • 10^°) . iHoranvi et al\ (|l999t) made mea- 




by Niemela and Sreenivasa 3 (|2003al ). iNikolaenko etali 
(|2005l ) where many data sets have been compared. There 



surements with liquid sodium (Pr = 0.005, Ra < 10®). 
Together with the results for helium gas, air (Pr=:0.7), 
and water (4 < Pr < 7) , these low-Pr data imply a strong 
increase of Nu with Pr at constant Ra, as shown in figure 
O top. For Pr larger than about one a saturation sets 
in and Nu becomes Pr-independent for some Pr-range. 
Recent results usin g helium gas at low temperatures 
(|Roche et al\ (|2002t) ) and covering the range 0.7 < Pr < 
21 suggest a very mild if any incre ase with Pr. Results ob - 
tained with vario us organic fluids (lAhlers and Xul (l200lD . 
IXia et al\ (|2002[ )) for Ra = 1.78 x 10^ and 1.78 x 10^ 
are shown in fig. [5l bottom, and indicate a maximum in 
Nu(Pr) near Pr ~ 3, followed by a very gradual decrease 
of Nu with Pr that can be described by Nu oc Pr"° °^ 
over the Pr-range of the experiments. 
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most of tlic features of Nu(Pr) observed in experiment. 
Wlien Ra is not too large, it predicts Nu ~ Pr^^* at con- 
stant Ra for Pr < 1, a maximum near Pr = 3, and the 
very gradual decline for larger Pr. For large Pr the GL 
prediction is shown by the solid lines in figure [5)3. Al- 
though the parameters of the model had been adjusted 
using data for Pr up to about 30 (including the open 
circles in the figure) , the model agrees with the measure- 
ments up to Pr ~ 2000. The large Pr behavior result- 
ing from the GL theo r y is d iscussed in more detail in 
iGrossmann and Lohsd ( 2001 ) and the small Pr behavior 



Grossmann and Lohse ( 20081 ) 



E. The aspect-ratio dependence of the Nusselt number 



Several experimeii ts (lAhlers and Xul (l200ll). 
iFleischer and Goldstein! (l2002h. iFunfschilling et al\ 



( 2005)^ iNiemela and SreenivasanI 



Nikolaenko et al. 



Wu and Libchaben 




iSun et al\ 



( 2006a ). 



2005aD, 



IXu et al\ (|2000f) ) have 



probed the dependence of Nu at constant Ra and Pr on 
r. Using water with Pr ~ 4, it is found for P < 5 that 
Nu increases, albeit only very slightly, with decreasing 
P. For larger F the measurements up to F = 20 suggest 
no further change, indicating that a large-F regime may 
have been reached. The weak F dependence suggests 
an insensitivity to the nature of the LSC (see also 
Sec. lin.F[) . which surely changes as F increases well 
beyond one, and is consistent with the determination of 
Nu by instabilities of the thermal BLs. - Theoretical 
efforts to understand the influence of F on Nu have been 
quite limited, see e.g. IGrossmann and Lohs3 ( 2003f ) and 
IChing and (120061 1 



F. The insensitivity of the Nusselt number to the LSC 

Several experiments suggest that the Nusselt num- 
ber in the Ra range below the transition to the ulti- 
mate regime is remarkabl y insensitive to th e strength 
and structure of the LSC. TCioni et al\ ( 19961 ) measured 
Nu(Ra) with a sample of water with Pr ~ 3 in a container 
of rectangular cross section in which the azimuthal LSC 
orientation was more or less fixed. They determined the 
heat flux both of the original water samples, and of the 
same samples after several vertically positioned screens 
had been installed within them. In the absence of the 
screens shadowgraph visualizations showed that plumes 
generated at the bottom boundary layer were swept later- 
ally just above the boundary layer by a LSC. The plumes 
rose vertically in the presence of the screens, suggesting 
a dramatically altered and much weaker LSC. For both 
cases the heat current was the same within a few percent. 
This experiment suggests that the heat current is deter- 
mined primarily by the conductance and instability of the 
thermal boundary layers which are not influenced signifi- 



enthalpy will find their way to the top on e way or another 
regardless of any LSC. ICioni et "oil (Il996h also found that 
tilting their cells relative to gravity by an angle a as large 
as 0.06 radian, which enhances the Reynolds number of 
the LSC, had no influence on the heat transport within 
their resolution of a few percent. 

More recentlv lAhlers et al. I (|2006b[ ) measured Nu(Ra) 
for a cylindrical water sample with F = 1 and Pr = 4.4 
as a function of the tilt angle a with a precision of 0.1%. 
They found, for example, a very small reduction, by 
about 0.4%, for a tilt angle a = 0.12 radian. In the same 
experiment the LSC Reynolds number was determined 
and found to increase by about 25% for Ra = 10^ and by 
about 12 % for Ra = 10^^. If the Reynolds number had 
any direct influence on Nu, one would have expected an 
increase of Nu with Re. Again one is led to conclude that 
the heat transport is independent of the vigor of the LSC 
and thus presumably determined by LSC-independent 
boundary layer properties. This finding seems to be in 
conflict with the flnal GL results eqs. (|2T|) and (|22|l . in 
which Nu and Re are intimately coupled to each other. 



For a F = 0.5 water sample IChilla et ai\ (|2004al ) 
measured a reduction of Nu by about 5% when they 
tilted their system by about 0.03 radian. Samples 
of this aspect ratio are more complex because the 
LSC can consist of either a single convection roll, 
or of a more complex structure approximated by two 
rolls s ta cked one ab o ve the other ( Verzicco and Camussil 
(l2003l ). IXi and Xial (|2008bf )V The authors conjecture 
that the tilt stabilizes the single-roll structure, and that 
this structure gives a smaller heat transport than the 
two-roll structure, thus accounting for the reduction of 
Nu. However, it seems surprising to us that for F = 0.5 
the Nusselt number should be more sensitive to the LSC 
than it is for the F = 1 system. 

More evidence for t he insensitivity of N u to changes 
in the LSC is given bv IXia and Luil ( 1997t ). who altered 
the LSC into an oscillating four-roll flow pattern by plac- 
ing staggered fingers on the side- wall and found that Nu 
changed very little. IXia and Qiul (|l999f) made an even 
stronger perturbation to the system by placing a baffle 
at the cell's mid-height, again finding insensitivity of Nu. 

In addition to the evidence of the insensitivity of 
Nu(Ra) to changes in the LSC, there is good evidence 
for the sensitivity of Nu(Ra) to the structure of the 
thermal BLs. This is pr ovided by an experiment of 
iDu and Ton3 (|2000l . l200lh who covered the top and 
bottom plates with triangular grooves that were much 
deeper than the BL thickness. They found an enhance- 
ment of Nu(Ra) by as much as 76%, with no signifi- 
cant change in the dependenc e on Ra (however, see also 
ICiliberto and Laroch i (I1999D 1. Flow visualization re- 
vealed an increase of plume shedding by the protrusions 
as the mechan ism of the Nu enhancement. Si milar results 
were found by it ringano and Verzicco! (120061 ) in their nu- 
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G. The dependence of Nu on Ra at very large Ra 

Below the transition to the uhimate regime the Nus- 
selt number is determined essentially by properties of the 
top and bottom thermal boundary layers (see Sec lIILFjl . 
As already discussed in subsections III.BI and III.C| this 
is expected to change dramatically in a critical range 
around some Ra* , defined by the condition that the shear 
across the laminar (albeit fluctuating) kinetic BL due to 
the LSC becomes so large that a transition to turbu- 
lence is induced within it. Note that the exact value 
of Ra* depends on the strength and type of the tur- 
bulent noise that perturbs the BLs, but the transition 
is expected to happen once the shear Reynolds num- 
ber Rcs, based on the kinetic BL thickness, exceeds 
Re* = 0(400). For r = 1 estimates of Ra * based on 
the GL theory (jGrossmann and Lohsd ( 20021 )) and cor- 
responding to Re* = 440 and 220 are shown in Fig. [Sb 
as dotted and dashed lines, respectively (since the pa- 
rameters of the GL theory have been determined only 
for F = 1, an equivalent prediction of Ra* for general 
r unfortunately is not available). These estimates are 
based on the assumption that a LSC continues to exist 
at these very large Rayleigh numbers. If it does not, 
then the transition should eventually be triggered by a 
destruction of the kinetic BL by turbulent fluctuations 
rather than by a laminar (albeit fluctuating) flow across 
the plates. Understanding the regime above Ra* is of 
particular importance because it is believed by many to 
be the asymptotic regime that permits, in principle, an 
extrapolation to arbitrarily large values of Ra, including 
those of astro- and geo-physical interest. 

Experimentally it should be possible to observe the 
predicted transition by a dramatic change in the mag- 
nitude and/or the Rayleigh- number dependence of the 
Nusselt number. For Nu(Ra) one expects a change from 
an effective power law with 7e// ~ 0.32 as observed be- 
low Ra* to 7e// ~ 0.4, which due to the logarithmic 
corrections is somewhat below the predicted asymptotic 
value 7 = 1/2 (see Sec. lILCj) . Another dramatic change, 
according to the theory, should be the dependence on Pr. 
For Ra < Ra* the Nusselt number is essentially indepen- 
dent of Pr for Pr > 1. For Ra > Ra* the Kraichnan 

prediction is Nu ~ Pr^^'''' (see Eq. [55]), at least for Pr 
near one. However the GL theory predicts Nu ~ Pr^''^ 
(see Eq. \T7} . so there remains some uncertainty on this 
issue. Nonetheless, any significant Pr dependence would 
lead to a discontinuity of Nu(Ra) of a size that would 
depend on Pr. 

From Fig. [HI bottom, one sees that the measurements 
with water at Pr = 4.4 and F = 1, 0.67, and 0.43 have 
not reached the regime above Ra* predicted for F = 1. 
As expected, the measurements for F = 0.67 and 0.43 
shown in Fig. [51 top, as well as those for F = 1 shown in 
Fig.m give no indication of the BL-turbulence transition. 
Neither do other measurements for Pr w 4, F between 



iNiemela et~al\ (|2000aD for F = 0.5 are shown in Fig. [3 
top, as solid circles. The data were corrected recently by 
some o f the original authors I Niemela and Sreenivasaiil 
(|2006a[ )) for side-wall and plate effects (see Sec. IlILB]).^ 
The data are for Rayleigh numbers as large as 10^^, 
and for Ra up to about 10^^ they are for Pr ~ 0.7. In 
the Ra-range of overlap, they are in excellent agree- 
ment with the water measurements for Pr = 4.4 and 
F = 0. 67 and 0.43, and with the resu lts for compressed 
gases ([Fleischer and Goldstein! (|2002[ )) with Pr ~ 0.7, 
1 < F < 3, and 1 X 10^ < Ra < 2 X 10^'^, demonstrating 
again the insensitivity of Nu to Pr and F in this Ra 
range, as well as a consistency between the cryogenic 
and room-tempera ture experiments. A t large Ra the 
values of Pr for the INiemela et al. I (l2000al) data increased 
because of the proximity to the critical point. As seen 
in Fig. [5)3, the F = 0.5 data might have been expected 
to cross Ra* somewhere near Ra = 10" or 10^'', but 
apparently did not do so since they reveal no change of 
the dependence of Nu on Ra. Two possible explanations 
come to mind. Perhaps the LSC was less vigorous in 
this experiment than it was for F = 1. In that case 
the expected Re* ~ 400 would only be reached at even 
higher Ra. Alternatively, at the very high Ra the LSC 
may have deteriorated into an unrecognizable entity 
consisting es sentially only of vigoro us fluctuations as 
suggested by ISreenivasan et al. 1 (120021 ). In that case the 
GL estimate for Ra* would no longer be quantitatively 
applicable. 

Within their experimental uncertainty and over the 
very wide range 10'' < Ra < 5 x 10^^ the data of 
INiemela et "oil (l2000aD can be described by a single power 
law Nu = A^oRa'^"^^ with an effective exponent jeff = 
0.323 and Nq = 0.0783. This power law agrees ve r y wel l 
with a fit to the data of [Fleischer and Goldstei 3 (I2OO2D 
which (over their much more narrow range of Ra) yielded 
Nq = 0.0714 and 7e// = 0.327. Both sets of measure- 
ments are inconsistent with an exponent of 1/3, which 
would yield a slope like that of the dotted line in Fig. [5^. 
However, they are remarkably consistent with the pre- 
diction of GL, which is shown by the solid line in the 
figure. The drop below the power law for Ra > 5 x 10^^ 
is unexplained. One might have attributed it to non- 
Boussinesq effects, but for gases near the critical point 
the se would cause an i ncrease of Nu and not a decrease, 
see [Ahlers et all (|2008l) and Sec. IVIII Alternatively, one 
might look at the variation of Pr as an explanation, but 
in the GL model Nu is essentially independent of Pr at 
these large Ra, cf. regime IVu in subsection III.BI 

An earlier set of data using helium at low tempera- 



0.67 and 20, and Ra up to 5 • lO^^ (|Sun et al\ (|2005al )) 



^ Interestin gly, this side - wall co rrection, which was based on the 
model of iRoche et al\ l l2001al ). yielded corrected data that are 
n early identical to those that had been obtained using Model 1 
of lAhler j l|200(]| 'l and had been shown in Fig. 5 of that reference. 
The plate-effect corrections are quite small for the cryogenic data 
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tures and F = 0.5 was obtained bv lChavanne et "oil (119971 
ll996Ll200lh. It extends up to R a ~ 10^^, and the results 
listed bv lChavanne et al\ ( 200l[ ) are shown as open circles 
in Fig. m In the range 10^° < Ra < 10^^ they agree very 
well with the other data sho wn in the figur e . For s maller 
Ra they are higher than the INiemela et al. I (l2000aD data. 
A possible reason might be found in a difference of the 
side-wall correction that was applied.^ More interesting 
is the difference between the two data sets that evolves 
as Ra grows beyond lO^'^. In that regime the open circles 
in the figure can be described within their scatter by a 
power law with 7e// — 0.38. The authors interpret this 
result as corresponding to the expected 7 = 1/2 in the 
Kraichnan regime, modified by the logarithmic correc- 
tions that are attributable to a viscous sublayer. Thus 
they claim to have entered th e "ultimate" , or asyrn p- 
totic, regime of turbulent RBC ( Chavanne et all ( 19971 )). 
However, the transition at Ra* just above 10^^ is lower 
than the theoretical estimates for the shear-flow bound- 
ary layer instability. For Pr = 1 and Ra = 3 x 10^^ the 
GL model yields (in the F = 1 case) Rcs ~ 100, which 
is too low for a shear-induced transition to turbulence 
in the boundary layer. An explanation in terms of a 
shear-induced BL transition would require a more vigor- 
ous LSC for the F = 0.5 case than was measured (see 
Sec. IIVI) for the F = 1 .0 case.* In any case, the data of 
IChavanne et al\ ( 200ll ). and the interpretation in terms 
of a transition to the Kraichnan regime, differ dramat- 
ically from the measurements of INiemela et 'aH (|2000aD 
who did not find this transition even though their data 
extend to higher values of Ra, and were done at some- 
what lower Pr where the shear transition should occur at 
even smaller Ra. The reason for this difference remains a 
mystery at this time, and the resolution of this apparent 
conflict between the two data sets is one of the major 
challenges in this field of research. 



Yet another set of data, shown as solid squares in 
Fig. E^, was obtained wi t h low- temperature helium by 
iNieniela and Sreenivasa 3 (l2003al) . using t he original ap- 
paratus of the F = 0.5 measurements by INiemela et al\ 
(|2000ar i. but with a sample of reduced height that had 
F = 1.0. Unfortunately these results do not help to clar- 
ify the situation. They are fairly consistent with other 



data in the Ra ran ge near 10 . At smaller Ra they agree 
fairly well with the lChavanne et~a B ('2001 ) data, but dif- 
fer from the side-wall-corrected iNiemela et al\ (|2000al ) 



3 In fact, for Ra < 10^ the INiemela etaB (|20003Ldata^.before 
their s ide- wall correction were much closer to the lChavanne et al\ 
(l200lh data than after this correction was made. 

* IChavanne et ai\ 1I2OOII ) measured the vertical LSC velocity com- 
ponent Vz{r,z,9) in the horizontal mid-plane (2 = 0) and at a 
radial position r = L/A half way between the center line and the 
side wall. The significant size found for Vz suggests the existence 
of well developed up-flow and/or down-flow in the mid-plane, im- 
plying that the LSC consisted of a single convection roll rather 
than of two rolls one above the other. The measurement yielded 
a Reynolds number that is not very different from that for the 
r = 1 case (sec Sec. ilVjl . but it is difficult to know precisely 
the corresponding shear across the viscous BLs at the top and 



data. Here again one would be tempted to invoke the 
side-wall effect as a possible explanation. More difficult 
to disregard are the data for Ra > 10^^, where side- wall 
corrections are negligible. There the data fall between 
the two F = 0.5 cryogenic data sets, thus adding to the 
complexity of experimental information about a possible 
Kraichnan transition. 

The Grenoble/Lyon group undertook several investi- 
gations in an attempt to find an explanation for the dif- 
ferences betw een the various data sets for Ra > 10^^. 
For instance, IChiUa et al\ (|2004bf ) developed a model 
that attempted to explain the difference in terms of a 
finite plate-conductivity effect (see Sec. IIII.B[) : but mea- 
su rements with relati vely low-conductivity brass plates 
bv lRoche et al\ ( 2005[ ) yielded results comparable to the 
high-con ductivity co pper-plate results. In a separate ex- 
periment iRoche et al . (2001b) made measurements using 
helium in a sample cell with F = 0.5 with walls and plates 
that were covered completely by grooves. The depth of 
the grooves was stated to be less than the thermal bound- 
ary layer thickness. Such a geometry is asserted to re- 
move the influence of the sublayer which is responsible 
for the logarithmic corrections. For Ra > 10^^ this ex- 
periment yielded an exponent quite close to 0.5, consis- 
tent with the expected Kraichnan va lue of 1/2. How- 
ever, iNiemel^^nT^Sreeni^Ka^ ( 2006a[ ) pointed out that 
the BL thickness decreases with increasing Ra and be- 
comes comparable to the groove depth in the Ra-range 
of the measur e ments . In such a case the measurements of 
IDu and Tond ( 2000[ ) using a sample with grooves in the 
plates that were deeper than the BL thickness indicate 
that the pre-factor of an effective power law describing 
Nu(Ra) increases by a s much as 76% for deep grooves . 
Thus, it is suggest ed bvlNiemela and Sreenivasa T(|2006aD 
that the results of lRoche et al\ ( 200 Ibl ) might possiblv be 
due to a cross-over between rough surfaces with a groove 
depth less than the BL thickness to a regime where the 
groove depth is larger than the BL thickness. More work 
seems needed to resolve this issue. 

An interesting e xperiment re l ated to the Kraichnan 
regime was don e bvlGibert et al! ("2006|) , taking up earlier 
experiments by iPerrier et ah (.2002 ) and experimentally 
realizing th e theoretically su g gested homogeneous RB 
turbulence (ICal zavarini et ail (|2006t ). iLohse and Toschil 
(120031) '!. iGibert et al. (2006.) used a vertical channel with 
wide entrance and exit sections that avoided the influence 
of the thermal BLs on Nu. They found relationships for 
Nu and Re (based on the velocity fluctuations) consistent 
with eqs. (|27|) and (|28|) when they re-defined Ra in terms 
of an intrinsic A-dependent length scale proportional 
to the ratio of temperature-fluctuation amplitudes and 
the vertical thermal gradient, instead of using a sample- 
geometry-dependent and A-ind ependent length. The 
same scaling was also obtained bv lCholemari and Arakeril 

Anr\f\ rl lnnr»ol\ 1 J„i j 1 1 j. 1 n 
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in a vertical pipe. The flow was driven by an unstable 
density difference across the ends of the pipe, created us- 
ing brine and distilled water. Away from either end, a 
fully developed region of turbulence existed with a linear 
density gradient. With a Rayleigh number based on the 
local density gradient relations consistent with eqs. ([27]) 
and (|28p were found experimentally. 

Also local measurements of the h eat flux ca n resul t 
in the 1/2 expon ent predicted by iKraichnanI ( 1962^ . 
IShang et aZ.I (200S ). based on their earlier measurements 
(jShang et al determined a time-averaged 

local Nusselt number at fixed positions x and given by 



Nu(a;) = {u3{x,t)e{x,t))J{KAL-'^) 



(30) 



in a r « 1 cylindrical cell filled with water, in the 
range 10^ < Ra < 2 • 10^". At half-height they found 
Nu(a;) = l.SRa"-^''^^" "^ when x was close to the side- 
wall where the heat was primarily transported by plumes. 



and Nu(a;) = 3.5 • 10"^Ra' 



0.49±0.03 



when X was close 



to the cell center where the plume density was much 
less^. Th ese two different local scaling l aws had been pre- 
dicted by[G rossmann and Lohsd (|2004l ) ; they correspond 
to the two independent scaling contributions eg_pi = 
KA^L~^Nuside and eg^^g = kA'^ L~^lSlucenter to the ther- 
mal dissipation rate eg ^ kA^L^^Nu (see equation (fT^ l. 
The theory predicts that for Ra < 10^^ the heat trans- 
port via the flow close to the side-wall dominates, but 
around Ra « 10^"' the heat flux through the center and 
thus the background fluctuations take over. Interestingly 
enough, when extrapolating the measured power laws for 
the local Nusselt number in the center and c l ose to the 
side- wall towards larger Ra, also lShang et al. I (I2008D ob- 
tain a crossover around Ra w 10^*. - The spatial inho- 
mogeneity of the local heat flux also has been confirmed 
by th e numerical simulations of IShishkina and Wagned 
(|2007D . 



IV. EXPERIMENTAL MEASUREMENTS OF THE 
REYNOLDS NUMBERS 

A. Reynolds numbers based on the large-scale convection 
roll 

The geometrical features and dynamics of the large- 
scale circulation (LSC) depend on the symmetry and 
aspect ratio of the sample. Here we focus on cylin- 
drical samples with F ~ 1 because for them the flow 
geometry is relatively simple, the dynamics is very 



^ Note that in earlier work bv lChine et all ||2004| ) a less steep scal- 
ing of Nu(a3) with Ra was derived at the center of the convection 
cell. One possible source of discrepancy is that the association of 
the length scale resulting from balancing buoyancy and viscous 
forces with the thermal boundary layer thickness done in that 



rich, and the experimental studies are most exten- 
sive. Then the flow occurs in a near-vertical plane 
and yields near-elliptical stream lines, with the long 
axis of the ellips e slightly tilted relative to gravit y 
(see, for instance. iQiu all (|2004l) . ISun et all (|2005c[) . 
IVerzicco and Camussil ( 1999f )). as illustrated in the time- 
average velocity map Fig. [1] bottom. Depending on Ra, 
two small counter-rotating vortices positioned near the 
corners close to the minor axis of the ellipse may be 
m ore or less prominent. I n the 2D numerical simulations 
of lSugivama et~al. I (Hooi) they can be very pronounced, 
depending on Ra. 

Often it is said that the large-scale convection roll is 
driven by the localized volumes of hot or cold fluid known 
as plumes that are emitted from the top and bottom ther- 
mal BLs as a result of a BL instability, and no doubt 
these plumes play an important role. But in addition 
some of the heat current is conducted across the BLs, 
warming or cooling the fluid adjacent to them. This in it- 
self, without the presence of the plume inhomogeneities, 
would drive the flow as it does closer to the onset of 
convection where there are no plumes. Thus presum- 
ably the LSC derives its existence from a combination 
of these two heat-transport mechanisms with the rela- 
tive importance of each depending on Pr and Ra. and 
indeed it seems difficult to separate one from the other. 
For Pr ~ 4 it is known that, away from the BLs, this 
heat current leads to a destabilizing time-averaged verti- 
cal gradient of the azimuthal average of the temperature 
which is strongest near th e side wall while the interi or 
is more nearly isothermal ( Brown and Ahler^ ( 2007bl )V 
Superimposed upon this azimuthally-averaged gradient 
is the warm up-welling and cold down-welling current of 
the LSC which, when time averaged, leads to a near- 
sinusoidal azimuthal temperature varia tion with period 
27r and amplitude 6 at the side wall (Brown and Ahlerd 
(l2007bD ). a nd a near- linea r tempe rature variation along 
a diameter (|Qiu and Tond (|2001bD ). 

Although the speed of the flow shown in Fig. [7^ varies 
considerably with position, one might expect that the 
large central roll can be described by a unique turnover 
time T. Using L as a relevant length scale, one can define 
a Reynolds number 



Re 



LSC 



2L2 
TV 



(31) 



We would not expect Re^"^*^ to describe all aspects of the 
flow field; for instance the small counter-rotating vortices 
might require a different Reynolds number which might 
even have a different dependence on Ra; but the main 
features, for instance those predicted by the GL model, 
might be related to Re^'^'^. Further, it has been sug- 
gested recently that the actual path length of the LSC 
circulation varies with Ra because the shape of the flow 
field changes; this feature would introduce an additional 
Ra dependence of Re^"^*" because the length scale, set 
equal to 2L in Eg. (IHTI). would ii o longer be a consta nt 

A-\Ti 1- 1 n ; I ilor\r\ol-l\ Ici 1 -vi-l Anr\r\rh\ 
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To our knowledge there are no direct measurements 
of T. However, T was inferred from local measure- 
ments of various velocity components (see Fig. [7^ for 
an example) and the assumption of a constant circula- 
tion path- length proportional to L. The ti me-averaged 
maximum vertical velocity co mponent Vmax ( Lam et al\ 
(|2002[ ).[q iu and Tonal (|2001a[ )) near the region between 
the viscous boundary layers and the bul k of the system 
(some times known as the "mixing zone" ( Castaing et al\ 
(I1989D ) gave Re^'""'" oc VmaxL/v. Alternatively, the slope 
7„ in the sample interior away from the mixing zone of 
the time-averaged horizontal component u of the LSC 
velocity as a function of the vertical position along the 
sample axis (see Fig [TK) was used to defi i ie and deter- 
mine Re" = "iuL^/y (|Qiu and Tongl (|2001all200l ) (open 
circles in Fig. [Th) . 

An interesting property of the LSC is a torsional 
azimuthal oscillation mod e with fre quency /n = a-'n/27r 
Funfschilling and AhlersI (|2004[ ). iFunfschilling" et al\ 



correlation functions between two temperatures mea- 
sured with probes s epara t ed vertically by a s mal l distance 
I (.Castaing et al. ( 19891) . IChavanne 6^0?] (IT997 . 20011 
iNiemela eii'aLl (| 2001(08 ano et aLl (|l989^ ■ lTakeshita etM 



2008|)) that can be used to define Re" = 2L2/(,/ 



Well before the spatial nature of this mode was 
known, its frequency was measured in numerous 
single-po int determinati o ns of t he temperature or th e 
velocitv (ICastaing a/] (Il989l). ICihberto et di\ (Il996l) 



ICioni aZ.I (Il997l). iHeslot e< a/.l (Il987l). iLam aZ.I 



(l2002D.lQiu et al\ ( 2004D.lQiu and Tond (l2001allbl. 120021) 
iQiu et al\ (l2000l) Takeshita et all (|l996l ). IXin and Xial 
(|l997t ). IXin et al\ (|l996f )). both of which have an oscil- 
latory contribution provided the probe is not located in 
the horizontal mid-plane of t he sample where th e ampli- 
tude of this mode vanishes ( FunfschiUing et all ( 20081 )). 
Some of the single-point measurements yielded results 
for Re'^ th at were equal to Re " within experimental 
resolution (|Qiu and Tond (|2002[ )). as seen in Fig. [7|d. 
The reason for this equality is not known at this time. 
Some other experimental investigations indicated that 
there is a distinct dif ference between th e Ra-dependence 
of Re"""" and Re" (|Lam et al\ (|2002l) ) and that these 
Ra dependences change with Pr. Clearly there remain 
some unresolved issues.^ 



B. Reynolds numbers based on plume motion 

An estimate of another Re was based on the motion of 
plumes. When a cold or warm plume passes a local tem- 
perature probe, then it produces a positive or negative 
deviation of the local temperature from the mean. The 
local time- averaged vertical plume speed Vpi = l/tpi was 
thus determined from the peak location tpi of time cross- 



A reliable evaluation of the experiments is complicated by the 
fact that experimentalists often slightly tilt their samples so as 
to obtain a dominant LSC orientation for their measurements. 
Measurements of Re are r emarkably sensitive to the sample align- 
ment relative to gravity llAhlers et al. and in a tilted 



( 19961 )'). The measurements yielded Re""' ex VpiL/v. 
A similar technique used thermometers mounted in the 
side wall at the horizontal midplane. A single temper- 
ature sensor yielded a time auto-correlation (AC) func- 
tion with a broad peak at a del ay time corre s pondin g 
to the plum e turn over time Tpi ( Ahlers et al\ (|2006bl ). 
iBrown et ^ (|2007t )). Similarly, two sensors on opposite 
sides of the sample yielded time cross-correlation (CC) 
functions with broad peaks corresponding to half the 
plume turnover time Tpijl. These measurements gave 
Reynolds numbers Re^' = {2L/Tpi){L/v). Quite remark- 
ably, measurements indicate over a wide parameter range 
that Re"""" = Re" = Re^' = Re""' within fairly small 
experimental errors. This can be interpreted to mean 
that the plume circulation is slaved to the LSC, or vice 
versa, and that all of these quantities (where they agree 
with each other) yield a reliable representation of Re^^'^ . 

The cross-correlation functions for two thermometers 
at the mid-plane but on opposite sides had extrema that 
were negative, indicting a correlation between the warm 
rising plumes on one side and the cold fallin g plumes on 
the ot her. This tends to support the ideas of lVillermauxl 
( 19951 ). who suggested that a plume impinging on the BL 
causes an instability and an associated emission of a new 
plume of the opposite type. However, the width of the 
cross-correlation function extremum, which was not very 
different from 7^/, indicates that this is not a periodic 
process as had been originally suggested. The periodic 
signal, when observed in some experiments, presumably 
is due to the torsional oscillation discussed above and in 

Sect. Ivml 

For modest Ra, say Ra < 3 x 10^ when 
Pr ~ 4, the GL prediction ( Grossmann and Lohsd 
(|2002D }_ is in good agree m ent with expe ri ment al re- 
sults (jBrown et all (|200(f ). iQiu and TonS (|2002[ )) for 
Rc"(Ra,Pr), Re"(Ra,Pr), and Ref'(Ra,Pr) (see Fig.©. 
In that parameter range all three measured quantities, 
so far as they have been determined, agree with each 
other. In that regime both the prediction and the exper- 
iment can be well described by a power law Re oc Ra'''^^^ 
with a very slightly Ra-dependent effective exponent 
7e// ~ 0.44 (there also are some notable differences 
between th e pred ictions and some other measurements 
(|Lam et all (|2002D )). However, for Ra > 3 x 10^ the mea- 
surements reveal a relatively sudden transition to a new 
state of the system, with a Re^' that is larger than Re" , 
which in turn is larger than the GL prediction. These 
results agree quite we U with recent experimental results 
for ReP' obtained by [S un and xial ((2005.1 in the range 
7 X 10^" < Ra < 2 X lO^^. If the GL model still cor- 
rectly predicts the Ra-dependence of Re^"^*"^ in this pa- 
rameter range, then the experimental results for Re^' and 
Re" suggest the existence of a new LSC state in which 
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other and where the twist osciUation of the LSC is no 
longer synchronous with the LSC turnover time (see also 
IXi et all (|2006. )). The observation that Re^' > Re^"^^ 
would then indicate that the plumes rise and fall more 
rapi dly than the background flow of the LSC. How- 
ever, iNiemela and Sr cenivasan ( 2003a[) and lSun and Xial 
( 20051 ) suggested that it is more likely that the LSC 
evolves into a more complex flow structure where its 
dynamics can no longer be described quantitatively by 
a uniquely-defined Reynolds number. It is unclear at 
present whether the difference between this state and the 
one at smaller Ra will be found in the geometry of the 
flow, in the nature of the viscous boundary layers that 
interact with it, or in the nature and frequency of plume 
shedding by the thermal boundary layers adjacent to the 
top and bottom plates. 

The dependence of Re'^ on Pr has been i nvestigated 
over a wide range of Pr by iLam et 'al. 1 (I2002D . The data 
for Ra < 3 X 10^ are shown as squares in Fig. [5]:. Their 
agreement with the GL prediction is not as good as one 
might wish. However, at Pr c± 4, they al s o diff er from 
the measurements of Re^' by iBrown et al\ { 2007t ) . More 
work over a wide Pr range seems desirable. 



V. NU(RA,PR) AND RE(RA,PR) IN DIRECT 
NUMERICAL SIMULATIONS 

Direct numerical simulations (DNS) of Rayleigh- 
Benard flow have several advantages in comparison to 
experiments: (i) Any local or global quantity can be 
"measured" without interfering with the probe or hav- 
ing restricted accessibility, (ii) The boundary conditions 
can be chosen as in the idealized RB case, i.e., with ex- 
actly zero heat flux through the side-walls and constant 
temperatures at the top and bottom plates. Therefore, 
a modeling of side-wall and plate corrections (see sec- 
tion [TlLB]) is not necessary, (iii) The material properties 
of the fluid can be chosen at will. E.g., one can study flu- 
ids which exactly obey the Boussinesq approximation, or 
fluids which show temperature dependence of one mate- 
rial property only. By doing simulations with such fluids, 
theoretically suggested mechanisms can be tested. 

However, obviously there are also major disadvantages 
of DNS, the main one being the limitations in the obtain- 
able sizes of Ra and Pr. For fixed Pr 1 the required 
CPU time increases roughly as Ra'^'^^ log Ra. The spatial 
resolution requirements for therma l convection have been 
defined bv lGrotzbachI (Il98lll983f ). In some simulations 
good small-scale resolution, which is particularly crucial 
in the BLs, is sacrified for higher Ra. 

In this DNS part of the review we again will focus on 
global properties such as Nu and Re and the dynamics 
of the wind, omitting DNS resul ts for the smal l s cales. 
The RB simu latio ns by DeLuca et al\ ( 199Cl[ ). 



to make efficient use of pseudo-spectral methods. Note 
that free-slip side- wall BCs lead to a different large-scale 
circulation (LSC) pattern than no-slip side-wall BCs. 
In particular, the counter rolls in the corners of the 
cell a re suppressed. With this method, iDeLuca et al\ 
(|1990D carried out simulations up to Ra = 1.6 • 10" 
for Pr — 1, finding an effective Nu vs Ra po wer-law 
expon ent consistent with 2/7 at the largest Ra. IWernd 
(I1993D extended these 2D simulations to Pr = 7. 

3D simulations in the early 90s - also employing spec- 
tral methods - could only achi eve relatively small Ra , e.g. 
Ra = 6.5 ■ 10^ in the wor k bv lSirovich et al\ (|l989f) and 
iBalachandar et al. lll989l) or Ra = 6.3 • 10^ in the work 
by Christie and Domaradzkil ( 1992f ): both are too small 
to make statements on exponents for power-law depen- 
dences on Ra, and also because there are still coherent 
flow structures. 

One of the first 3D RB simulations fro m wh i ch suc h ex- 
ponents could be extracted was that bv lKeril (|l996[) . He 
employed lateral periodic BCs and achieved Ra = 2 • 10^ 
on a 288 x 288 x 96 grid. For Pr=0. 7 he found an ef- 
fective power law Nu ~ Ra^ '^^. Later, iKerr and Herring 
(I2OOOD extended these simulations to Prandtl numbers 

for 



Ra' 



0.25 



in the range 0.07 < Pr < 7, finding Nu 
the smallest Pr and effective exponents consistent with 
2/7 for Pr > 0.7. For the Prandtl dependence of Nu 
at fixed Ra = 1 ^ and Pr < 0.7 (see figure [5l top), 
iKerr and Herrind (|2000[ ) give an effective exponent of 
0.12. A m ore recent example o f a spectral RB code is 
the work of lHartlep et al. achieving Ra = 10''. In 

that work the focus was on the flow organization in RB 
convection with sidewise periodic boundary conditions. 

A second class of DNS for RB convec tion is that 
of La ttice-Boltzmann (LB) simulations ( Benzi et all 
(I1994D1. In one o f the first large LB RB simulations 
iBenzi et all (|1998D achieved Ra « 3.5 • 10^ and found 



7iVii = 0.283 ± 0.003 for the effective Nu vs Ra power- 
law exponent (they used Pr = 1 and free slip on vertical 
walls). 

The third class of simulations contains DNS 
based on finite-difference or finite-volume. 

their 



IVerzicco and Camussi 


(I1999D. building 


earlier work 


Camussi and Verziccol 


IVerzicco and Camussi 


(19971)). emploved 



on 



(IT999I) . 



IWerne et ail (|l99l[ ). and lWernd (l993[) were restricted to 



two dimensions (2D) and employed free-slip (but imper- 



difference simulations to obtain effective power-law 
exponents for both Nu and Re vs Ra for several Pr for 
Ra up to 2 • 10^ (see figure [5^). The advantage of this 
method is that the correct no-slip BCs at the side-wall 
can easily be implemented and that simulations can also 
be performed efficiently in the cylindrical geometry as 
used in many experiments, so that a direct comparison 
with the experimental results is possible. A further 
advantage of finite-difference or finite-volume methods 
is that the spatial resolution can easily be refined at will 
in the BLs, using a coarser grid in the biflk. For the 
small Prandtl number Pr = 0.025 and F = 1 the effective 
power law Nu w 0.1 WRa*^'^^ wa s found in the regime 
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(|l999l ) obtained the effective power law Nu - Ra°-^^^. 
These authors also explored the Pr dependence of Nu 
over the range 0.0022 < Pr < 15 at fixed Ra = 6 • 10^, 
finding Nu w 8.5Pr°'^'' as an effective power law. For 
the Reynolds number they found Re ~ pj.-0-73 
Pr < 1 and Re ~ Pr-"'^'' for Pr > 1. Similarly, 
iBreuer et al\ i 20041 ) performed finite volume integrations 
of the Boussinesq equations, finding Nu ~ p^o.i82±o.oi2 
and Re - p^-o.607±o.oi2 f^^. ^q-s < Pr < 1 and 

Nu - pr0.032±0.003 aj^d _ p^-0.998±0.014 fo^ 

1 < Pr < 10^, for no-slip boundary conditions and 
R a = 10^. 

Verzicco and Camussil ( 20031) and 

IStringano and Verzicc d (I2006D extended the earlier 
simulations to the remarkably high Ra = 2 • 10^^ for a 
slender F = 1/2 ceU and Pr = 0.7. Beyond Ra = 10^ 
the numerical data are consistent with Nu ^ Ra}/^. 
The focus of those papers is on the flow organization 
in the F = 1/2 cell: Beyond 10^° and for Pr « 0.7 the 
single convection roll can break up into two smaller 
counter-rota ting rolls, each appro x imate ly of aspect 
ratio one. IStringano and Verziccol ( 2006f ) in addition 
showed that the thermal properties of the side- walls can 
stabilize the large-scale convection roll. 

At present th e large st-Ra D NS of RB flow are the ones 
by Amati et al\ ( 2005f ) and bv lVerzicco and Sreenivasaiil 
(l2008|), achieving Ra = 2 • lO" for F = 1/2 and Pr = 0.7. 
The focus here is on the difference between constant- 
temperature and constant-flux boundary conditions at 
the top and bottom plates, which already have been dis- 
cussed in subsection IIFBI Similarly, side-wa ll correc- 
tions (|Verziccol (2002)) and plate corrections (jVerziccol 



( 2004 )') have been studied, which also already have been 
discussed in the same above subsection. 

Another advantage of finite-difference simulations is 
that complicated geometries like those w ith rough walls 
can be treated, see IStringano et all ( 2006[ ). The results - 
an enhanced heat flux consistent with Nu ~ Ra°'^^ over 
the range 2 • 10^ < Ra < 2-10^^ - are in reasonabl e agree - 
ment with the experimental results of lQiu et al. I (Hooi), 
who found a power law exponent of 0.35 in the range 
10^ < Ra < 10"'^°, but in conflict with the earlier men- 



tioned results bv lDu and Tond ( 2000| ). who found an un- 
changed effective power law exponents f or the rou g h wal l 
case, but a larger prefactor. Similarly, IVerzicccI ( 20031 ) 
found that by manipulating the velocity BC at the plates, 
the viscous BL could be affected and also the absolute 
value of Nu, but its power-law exponent was rather ro- 
bust against such manipulations. Clearly, more research 
is necessary to clarify this matter. 

The numerical s imulations of IShishkina and Wagneil 
(|2006l . l2007l l2008f) also are based on a finite volume 
scheme; they focus on the role of plumes and the flow 
organization in RB flow with Pr = 0.7 and Pr = 5.4 and 
Ra up to 2 • 10^. For Ra = 2 • 10^° large-eddy simu- 
lations (LES) were done. In this subgrid modeling the 
length scales in the dissipative and diffusive regime are 



Finally, when full spatial resolution of the turbulence 
field was abandoned even further, obviously much larger 
Ra could be achieved. An example for such a calcula- 
tion is the so-called transient Rey nolds averaged Navier- 
Stokes (RANS) method, which iKenieres and Hanialid 
(j2002i ) applied to RB convection. These authors find 
Nu - Ra"-^^ between Ra = 10^ and lO^^. For the last 
two Ra decades some reorganization of the plumes and 
a slightly enhanced Nusselt number were observed, but 
given the progressively decreasing spatial resolution of 
the numerical scheme at these high Ra the implications 
of this finding are presently unclear and a detailed discus- 
sion of RANS simulations and LES of RB flow is beyond 
the scope of this review. 

Given the enormous CPU power needed to achieve 
large Ra, one may wonder whether two-dimensional sim- 
ulations would not be sufficient to reflect at least some 
aspects of the dynamics of the three-dimensional RB 
problem. Thi s point has b een analyzed in detail by 
ISchmalzl et al. I (I2OO2I |2004^ whose conclusion is that 
for Pr > 1 various properties observed in numerical 
3D convection (and thus also in experiment) are in- 
deed well reflected in 2D simulations. This in particu- 
lar holds f or the BL pr o files a n d for the Nussel t num- 
ber. Also lAhlers et aZi (|2008[ ). ISugivama eFall (|2007l 
I2OO8I ) employed 2D numerics to study the non-Oberbeck- 
Boussinesq deviations of the Nusselt number and of the 
bulk (central) temperature from that in the Oberbeck- 
Boussinesq cases, see section IVlIl 



VI. BOUNDARY LAYERS 

A. Relevance of boundary layers and challenges 

Boundary layers describe the temperature and flow 
fields in the vicinity of the plates and walls. They are 
characterized by their time-averaged profiles in the di- 
rection perpendicular to the respective solid boundary, 
i.e., in the z-direction off the bottom and top plates or 
in the x-direction off the side walls. As the top and 
bottom boundary layers contribute the main resistance 
for the heat transfer through the cell and thus domi- 
nantly determine the Nusselt number, they deserve spe- 
cial attention. Indeed, nearly all theories of Nu(Ra,Pr) 
in RB convection are in essence boundary-layer theo- 
rie s: This holds fo r the now classical mixing-layer theory 
of ICastaing et al. the turbulent-BL based scal- 

ing theory of Shraiman and Siggial ( 1990f ) and extensions 
thereo f (e.g. Chiii^ 1997)), the turb ulent BL type theo- 
ries of lDubrullel (|200ll . |2002[) and of iHolling and Herwig 
(12005, 2006|)^ 

and the theory of iGrossmann and Lohsd 
(12000. . , 2OOIII2OO2II2OO4D which scaling- wise builds heavily 
on the laminar, though time dependent Prandtl-Blasius 
BL theory as elaborated in subsection III.BI It is there- 
fore of prime importance to directly study the BLs, to 
see whether the assumptions on which these theories are 
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Unfortunately, directly studying the BLs in the large- 
Ra regime is equally challenging experimentally, numer- 
ically, and theoretically. Experimentally, the required 
spatial resolution is very difficult to achieve. Even for 
the L « 6.5 m high Ilmenau RB barrel both the thermal 
and the kinetic BL thickne sses are only a f e w milli meters 
when e.g. Ra « 10^^^ see Idu Puits et al\ (|2007al [bh. In 
addition, both laser-Doppler velocimetry (LDV) and hot- 
wire anemometry are notoriously difficult to employ in 
regions with large temperature fluctuations and a small 
mean- flow velocity; both complications occur in the BLs. 
Numerically, large Rayleigh numbers are hard to reach, 
in particular with sufficient spatial resolution, see sec- 
tion IVl Even state-of-the-art DNS such as those of 
IVerzicco and SreenivasanI (2Q0l) have only a few grid- 
points in the BLs. Theoretically, no generally accepted 
BL theory for flow over a strongly heated surface exists. 
Even in the truely laminar case an analytic theory for 
the two-way coupling of the temperature to the velocity 
field is missing (i. e., taking the temperatu r e as an ac- 
tive scalar), see e.g. ISchhchting and GerstenI (|2000t ). And 
even if the two-way coupling is suppressed as done i n the 
interesting nume r ical w ork by IChing and Lol ( 200l[ ) and 
IChing and Pand ( 2002| ). no exact analytical results ex- 
ist for the thermal BL thickness and thus for the heat 
transfer through the shear flow over the heated plate. 

The analysis of the BLs is further complicated by their 
extreme complexity, not only in time, where there is 
plume detachment (see figure [1]), but also i n space. Re- 
cent ex perimental studies such a s those by Lui and Xial 
(119981). iMavstrenko e^ al\ (|2006D . IWang and Xial l{20m . 



and 



Sun et al. ( 20081) and numerical studies such as those 



bv ISugivama et al\ (|2008f) give increasing evidence that 
the BL thicknesses as well as the profiles of the mean 
quantities and of the fluctuations depend on the posi- 
tion along and above the plate(s), not only relative to 
the walls, but also relative to the fluctuating large con- 
vection roll (see section ITV)) . It is this main convection 
roll which creates the viscous BL because of the no-slip 
condition at the top and bottom plates. 

The aim of this section is to give an overview of what 
has been found experimentally, theoretically, and numer- 
ically about the BL thicknesses. 



B. Thermal boundary layers 

The thermal boundary layer thickness can be defined 
in several ways. From an experimental point of view, 
it is easiest to time-average the temperature profile at 
fixed lateral positions (cc, y) = * for various z and to ex- 
tract a thermal boundary layer thickness Xg{*) from the 
resulting local profile. From the theoretical viewpoint 
area-averages on x, y (in addition to the time-averaging) 
are of main interest. Indeed, it is the thickness \g of 
the area- averaged temperature profile in the z-direction 
which enters theory while comparison with experiment 



Note that strictly speaking the global thickness Xg is not 
equal to (Ae(*))^ in general; the reason is that the thick- 
ness is a (nonlinearly defined) property of a given profile. 
Then Xg in particular is a property of the area-averaged 
profile and not the average of properties Xg{*) of the re- 
spective local profiles. Along large plates the local BL 
width may well represent the area-averaged global one; 
but for RB samples with F of order 1 one has to expect 
significant differences, even different scaling of the global 
and the local thicknesses, since the area-averaged profiles 
comprise all inhomogeneities of the temperature and flow 
fields along the plates, including backflows, the regions 
between them, and the near wall ranges. These all will 
locally have quite different profiles. The more homoge- 
neous the flow, the more similiar the global width will be 
to the average of the local widths. Experimental analysis 
of all these details is still a big challenge. 

The (global) thermal BL thickness Xg can be defined 
in several ways. The most popular way is to define 
Xg through the slope of the area-averaged time-mean 
temperature profile at the plates: Take the tangent of 
the area-averaged mean temperature 9{z) at the plate. 
That distance between the plate and the vertical posi- 
tion where this tangent crosses the bulk (or center) tem- 
perature Tc is then called Xf . The center temperature 
Tc is equal to the mean temperature Tm ~ (Tt + Th)/2 if 
the Oberbeck-Boussinesq approximation holds because of 
up-down-symmetry. The very notion of Tc already refers 
to an area-averaged profile. 

Locally, i.e., for fixed {x,y) = *, the slope of the time- 
averaged profile, as well as the temperature for z oo 
and thus the local width Xf{*), will depend on the hori- 
zontal position * where the local time-mean temperature 
profile is taken. 

One of the first measurements of temper ature and also 
veloc ity profiles in RB cells was done by iTilgner et al\ 
(|l99t^l . namely in water (Pr = 6.6) at fixed Ra = 
1.1 • 10 ^ and at a fixed lateral position *. iBelmonte et all 
(I1993D extended these measurements to the Ra-range 
5 • 10^ < Ra < 10^^ in compressed gas (air) at room 
tempe rature (Pr = 0.7), but still at fixed lateral posi- 
tion. iLui and Xial ( 19981 ) measured the lateral depen- 
dence of Xf{*) on the positions x (in the mean LSC di- 
rection) and y (perpendicular to the preferred LSC di- 
rection) in a cylindrical water filled RB cell with F = 1 
in the regime Ra = 2 • 10® — 2 • 10^°. They report vari- 
ations of the BL width by nearly a factor of 2, depend- 
ing on the lateral loca t ion w here the profile was mea- 
sured. IWang and Xial (|2003') found similar results for 
a cubic cell. Even the scaling exponent of Xf{*) with 
Ra depends on the posit i on an d varies between —0.35 
and -0.28. iLui and Xial (|l998l ) found the thermal BL 
to be thinnest close to the center of the plates; there 
2Xf{0,0)/L = (0.23 ± 0.02)Ra-°-285±0-004 ^^^^^ 

One cannot draw conclusions about the dependence of 
Nu on Ra from such measurements of local Xf{*) at a 
particular lateral location. It is the length scale resulting 
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is connected with the Nusselt number through the exact 
relation 



Nu 



\dz {0{x, y, z = or z = i))^ | 



2\f 



(32) 



which follows from equation ([S]). 

A different definition for a thermal boundary-layer 
thickness is the one employed by cl assical ( la terall y 



1908), 



19611 ). 



homogeneous) laminar BL th e ories IIBlasiua 
Prandtll (ll905^■ see also l CowlevI (|200lh . iMeksvnl 
Schlichting and GerstenI ( 2000t )). This thickness, known 
as }?^'^° , is defined as the vertical distance from the plate 
to the point where 99% of the temperature difference be- 
tween plate and mean center temperature is achieved. 
E.g. for the bottom plate Ag^^° is the distance where the 
temperature reaches T = — 0.99Aft. This definition is 
in analogy to the common definition of the thickness 5 of 
the kinetic BL, defined by the distance where, say, 99% 
of the maximum bulk velocity is achieved. The profile- 
based thickness Ag^'^° can be calculated within the truly 
laminar Prandtl-Blasius BL theory, which implies lat- 
eral homogeneity and thus yields local and area-averaged 
widths that arc the same. One obtains 



A 



99% 



C(Pr) 



l/2p,l/3 



Re'/^Pr 



(33) 



with a hmction C(Pr) given bv lMeksvnl(|l96l[ ). For large 
Prandtl numbers C(Pr) becomes constant, whereas for 
small Pr one finds C(Pr) oc Pr"^/^. We note that Xf^' 
can display different scaling behavior than \f , both as a 
function of Ra and of Pr, depending on the parameter- 
space regime. To our knowledge this issue has not yet 
been explored systematically, neither for A^^'^" , nor for its 
local counterpart Ag^'^''(*) which depends on the lateral 
position * and refers to the local, only time-averaged, 
profiles. 

A thir d way to define a ther mal BL thickness, sug- 
gested bv lBelmonte et all ( 19941 ). is to take the position 
of the maximal temperature fluctuations at the edge of 
the thermal BL. We call the respe ctive BL thickn e sses de - 
fined in this way A^ and A^(*). Idu Puits etai\ (|2007b( ) 
show that thermal BL thicknesses defined in this way and 
measured in the Barrel of Ilmenau using air at one atmo- 
sphere as the fluid, behave differently from A^'. However, 
since BL thicknesses are determined largely by diffusive 
processes, it cannot be ruled out that a mixture like air 
(where both mass diffusion and heat diffusion play a role 
and where marginal stability depends both on concen- 
tration and temperature gradients) might in this respect 
behave differently from a pure fluid. 

Profile measurements, and thus the determination of 
BL thicknesses, are obviously less difficult if the abso- 
lute length scale L is larger. That is why the about 6.5 
m high ba rrel of Ilmenau i s partic ularly suited for these 
purposes. Idu Puits et al. I (l2007bD measured local time- 
averaged temperature profiles up to Ra ~ 7.7 ■ 10^^, for 



near the plate they obtain a power law 6{z) ~ z" with 
a w 0.5, only weakly dependent on Ra and F. For even 
larger z the temperature saturated at its bulk value. The 
1/2 power law for intermediate heights supports neither 
the logarithmic scaling laws for turbulent BL theories nor 
the temperature profile for a one-way coupled tempera- 
ture field (i.e.; taking the temperature as a passive scalar) 
predicted by the (laminar) Prandtl-Blasius approach. 
However, the power la w Xf{*)/L cx Ra~^^^^'^ found by 
Idu Puits et all ( "2007bl ) below the center of the top plate 
is not inconsistent with what follows from the scaling - 
wise laminar Prandtl-Blasius approach. ISun et all ( 20081 ) 
obtained Xf{*)/L oc jJa-0-33±o.05 fj-Qjj^ temperature- 
profile measurements above the center of the plate in 
water with Ra up to 2 • 10^°. 

Thermal BL thicknesses also have been determined 
from DNS. For a F = 1 RB cell IVerzicco and Camussil 
(I1999D obtained for the area-averaged width Xl = 
3.1i?a-0-29 for Pr = 0.7 over the range 5 • 10^ < Ra < 
2 • 10^ and A^' = 2.SRa-° '^^ for Pr ^ 0.022 in the regime 
5 • 10^ < Ra < 10^. When goi ng towards the larger 
Ra reg ime 2 • 10^ < Ra < 2 • 10^^ IVerzicco and Camussil 
([2003^ find XI - i?a-"-3i for Pr = 0.7 and F = 1/2. For 
even l arger Ra up to 2 • lO^'' IVerzicco and SreenivasanI 
(|2008l ) report A^ - Ra-^^^. 



The profiles of the rm s temperature fluctuations ar e 
much less settled, see e.g. Fernandes and AdrianI (120021) 



iLui and Xial (Il998l). Idu Puits et al\ (l2007bl). ISun et all 
(|2008l ). IWang and Xial (|2003l ). This issue will not be dis- 
cussed in the present review. 

As stated already in section III.B|, thermal plumes 
can be viewed as detached pieces of the thermal BL 
Since their discovery in RB convection by IZocchi et all 
( 1990l ) they have iDcen studied so intensively, - ex- 
perimentally, numerically, and theoretically ~ that a 
review of its own on this subject would be justified. 
A number o f recent papers were devoted t o them 



(iBreuer fflZ] (I2004D. i FunfschiUing and Ahlersl (I2004D 




iHaramina and Tilmed (^2004^. iPuthenveettil et all 



.Puthenveettil and Arakeril 
Shishki na and Wagned (12006 



Theerthan and Arakeril (Il998l. l20qol). Ixi et al 




Zhou et all (|2007bD . IZhou and Xial (I2002D). Their na 



ture w as discussed most recently by Funfschilling et "oZi 
(I2008D . Here we only make a few remarks. They 



seem to originate as one-dimensional excitations of 
the marginally stable BL. As they are born out of 
the BL, they become oriented by the LSC with their 
long axis in the direction of the flow, and then are 
swept by the LSC toward the side wall. Along their 
way they separate from the BL and progress vertically 
near the wall, developing th eir fa. mo us mushr o om to p 
(IPuthenveettil and Arakerl (l2005i). IXi et all (|2004l ). 
IZhou et all (|2007b[ ). IZocchi et all (|1990D ) in this process. 
Although they are often referred to as 'sheet-like", 
implying an extension in two spatial dimensions, the 
term "line-like" would more appropriately describe their 
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the thermal boundary-layer thickness in the other two. 



C. Kinetic boundary layers 

The kinetic BL thicknesses A„ are determined from 
the velocity profiles or the velocity- fluctuation profiles. 
Also A„ can be defined in various ways, starting either 
from time- and area-averaged quantities or only locally 
time-averaged ones. An additional alternative here is to 
average the velocity vectors themselves or to take the 
rms velocities, component wise or the full magnitude. 
Given the velocity profile of interest, its thickness may 
be defined (i) via the slope of the velocity profile in the 
z-direction (for the BL above or below the bottom or 
top plate) at a given position * on the plate, known as 
A^'(*), (ii) through the distance A™(*) to the local maxi- 
mal mean velocity, (iii) through the distance X^i*) to the 
maximal velocity rms fluctuations, etc. Again, all these 
quantities depend on the lateral position * = {x,y), un- 
less area-averages are considered. 

The first direct systematic measurements of velocity 
profiles in RB convectio n as a function o f Ra w here 
done by IXin et al\ (|l996f ) and IXin and Xial (|l997[ ). us- 
ing a novel light-scattering technique in a cylindrical cell. 
They found X^{*)/L ^ i?a^"'^^ from the velocity profile 
above the center of the lower plate. At the side-walls 
at half height as functions o f x instead o f z they mea- 
sured A"„'(*)/i - Ra"^-^^ . iQiu and Xial ((199813) ex- 
tended these measurements to convection in cubic cells, 
finding the same scaling exponents as in t he cylindri- 
cal ca se. By using various organic liquids, iLam et oil 
( 20021 ) could also explore the Pr dependence over the 
range 6 < Pr < 1027 and 2 • 10^ < Ra < 2 • 10^°, finding 



Pr°-24Ra- 



0.16 



above the center plate posi- 



tion. The small Ra-exponent of the effective power law 
is remarkable, as for a Prandtl-Blasius area-averaged BL 
profile one would expect an exponent of —0.25, just as 
measured at the side- walls, roughly corresponding to the 
Prandtl-Blasius scaling A^'/L ^ Re^^^^ with (approxi- 
mately) Re - Ra^/^ We will come back to this point 
later. Note that all these measurements were still done 
above a single position over the plate. 

In 2003, particle image velocimctry (PIV) measure- 
ments revolutionarized the experim ental analysis of th e 
veloc i ty field in therma l convection ( Sun et al\ (j2005bl l3. 
l2006l). (xia et al\ (l2003l)), inc luding the analysis of the ki- 



netic BLs (Sun et 



(200i)). This allowed not only for 
the direct identification of various Reynolds numbers, cf. 
Sect. IIVI but also for that of various kinetic BL thick- 
nesses such as A^'(*), A™(*), or A5^(*), even at different 
positions * in the cell. 

For 10^ < Ra < 10^° and Pr = 4.3 in a rectangular cell 
with r = 1, ISun et al\ jlooi) obtained time -averaged 
velocity profiles in the center above the bottom plate. 
They found that if the vertical lengths are rescaled by 
-^M (*) ^i^d the velocities by the maximal velocity, then 



on Ra; in this sense it is unive r sal iii that regime, see 
figure Eb. Moreover, [Sun_eLdJ 120081) find X'J{*)/L ~ 
K{*)/L ~ Ra-'^-^^=^°-°^ (see figure [Hi), in contrast 
to the much weaker R a dependence repo rt ed earlier 
(!_Qiu_ and xial (Il998allbl ). IXhT and xial 11993), IXin et al\ 
(199g). The origin of this discrepancy is not clear at 
this point. However, we note that the newly found scal- 
ing X'^{*)/L ~ A™(*)/L ~ Ra-^-^^^^''-^^ is consistent 
with the Prandtl -Blasius expectation. As a double check, 
ISun et al\ (|2008l ) plotted A^'(*)/L and A™(*)/i against 
the independently determined Reynolds number, find- 
ing the Prandtl-Blasius scahng Aj;'(*)/L ^ A™(*)/L ~ 
jjg-o.50±o.03^ see figure [TUb . In nu merical simu l ations 
of RB flow for Pr = 1 and F = 1 .Breuer et al\ (|2004l ) 



find X'^/L ~ Re' 



0.44±0.02 



, but a we aker Re-dependence for 



A™/L . Numerical simulations bv IVerzicco and Camussil 
( 19991 ) for the same F = 1 and Pr = 0.7 gave an area- 
averaged profile thickness exponent consistent with the 
Prandtl-Blasius BL th eory, namely Af.'/L = . 95j^a~ °-^'^. 
For lower Pr = 0.022 FVerzicco and Camussl ( 19991 ) ob- 
tained XfjL = 0.1i?a"0i8. 

The PIV BL study of ISun et al\ (|2008D also yielded 
the wall shear stress Tu,, the viscous sublayer length- 
scale j/iu = i'{z = 0)/ur, the skin- friction velocity 
Ut = (tw/p)^^^, and the skin- friction coefficient c/ = 
Tw/{pU^^^). (All these quantities have been measured 
above the plate center, but here for simplicity we suppress 
the * = {x, y)-dependencc in the notation.) They found 

Re°-^°, and c/ Re'^-^'^. The respective Prandtl- 
Blasius scaling exponents for a laminar BL over a flat 
plate are 3/2, —1, 3/4, and —1/2. Apart from the last 
one all are in fair agreement with the measurements. It is 
hard to judge the importance of the discrepancy between 
—0.34 and —0.50 for the scaling of the skin-friction co- 
efficient, in particular because the last three data points 
around Re « 10^ s how a stee p er dep endence than —0.34, 
see figure 22d of ISun et all (|2008|). Moreover in the 
numerical simulations of Iverzicco and Camuss I (I2003D 
(F = 1/2, Ra up to 2 • 10"'^^) the friction coefficient c/ 
does show a Reynolds-number scaling-exponent consis- 
tent with —1/2. (That work did not find any scaling 
for A^ as function of Ra, presumably because of a re- 
organization of the flow in the F = 1/2 cell with in- 
creasing Ra.) Also in the h i thert o largest numerical RB 
simulations of lAmati et all ( 20051 ) c/ ~ Re~^^^ is found, 
see figure 6b of that paper. Next, a lso the experimen- 
tal analysis of lChavanne et 'dl (120011) finds a —1/2 scal- 
ing of a surrogate of the friction factor with Re up to 
Ra = 10^^. However, at Ra = 10^^, corre s pondi ng to 
Re = 8 • 10*, the data of IChavanne et al\ (|200lD give 
evidence for a t ransition towards a weaker dependence. 
IVerzicco and Ca mussi (2003) find the same transition in 
this friction factor surrogate, but not in the actual fric- 
tion factor Cf. 

Finally, we wo uld like t o ment ion the interesting nu- 
merical study bv lYu et al. I (I2007D . who analyse the heat 
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steady convection with Ra up to 10^, finding basically 
the same Ra-scaling laws for Nu, A„, and Xe as in clas- 
sical RB convection. Their conclusion is that turbulence 
does not play a decisive role for the heat transfer. 

In summary, it is fair to state that, although some of 
the various velocity BL profiles differ from the idealized 
Prandtl-Blasius truly laminar profile due to the perma- 
nently ongoing plume emission, the scaling of the kinetic 
BL thickness with Ra, Re, and Pr is consistent with the 
laminar Prandtl-Blasius theory at least up to Ra = 10^^, 
but presumably also beyond. There is no indication of 
any transition towards a different kinetic BL thickness- 
scaling even at the largest Ra realized up to now. More- 
over, for Pr > 1 the thermal BL is nested in the kinetic 
one, whereas for Pr < 1 it is the other way round. 

The scaling of the r ms velocity fluctuations (see e.g. 
Idu Puits etdl (|2007al ). [Smi al\ ^OS) for recent work 
on this issue) will not be discussed in this review. 



VII. NON-OBERBECK-BOUSSINESQ EFFECTS 

The problem of RB convection is commonly analyzed 
within t he so-called Oberb e ck-Boussinesq (O B) approx- 
imation (|Boussinesd (|l903[ ). lOberbeckl l|l879t) ). in which 
the fluid properties are assumed to be temperature in- 
dependent, apart from the density for which the lin- 
ear temperature dependence eq. ([T|) is assumed. Un- 
der normal conditions, i.e., "small" temperature differ- 
ences A between the bottom and top plates, this ap- 
proximation is rather good. However, in order to achieve 
ever larger values of Ra for given L and fluid properties 
/?, Cp, K, and v, the temperature difference A between 



top and bottom plate quite frequently was increased to 
such an extent that the OB approximation had to be 
expected to fail. Non-Oberbeck-Boussinesq (NOB) ef- 
fects, i.e., deviations of various properties including Nu 
and the center temperature Tc from the OB case, then 
had to be expected at the larges t Ra i n several experi- 
ments ( Ashkenazi and Steinberd (Il999l). Castaing et al\ 



(Il989l). IChavanne et all (Il997l). iNiemela et all ll2000al). 



2006bD). Particularly 



iNiemela and SreenivasanI ( 2003 
problematic were measurements in the vicinity of liquid- 
gas critical points, where Ra tends to become exception- 
ally large but where NOB effects are to be expected at 
relatively small A. To be able to interpret these high- 
Ra data - in particular with respect to the question of 
whether there is an intrinsic large-Ra transition in the 
flow - it is therefore of prime importance to understand 
the physical nature, the signatures, and the size of the 
NOB effects. In addition to this practical consideration, 
there is much interesting physics to be learned from their 
study. 

NOB effects in high-Ra convection were measured and 
analyz ed first using helium gas by IWu and Libchabei 



a symmetry breaking between the top and the bottom 
of the sample. The temperature drops A^ and At across 
the bottom and top boundary layers become different, 
Afc 7^ At, and so do the thicknesses of the thermal BLs, 
^fb 7^ ^e\- Both phenomena are associated with a corre- 
sponding shift of the temperature Tc = Tf, — A^ = Tt -I- At 
in the center (bulk) of the cell away from the arithmetic 
mean temperature Tm = (?& -I- Tt)/2 of the bottom and 
top plate temperatures Tt and Tf Moreover, one would 
expect NOB observables to include deviations of Nu and 
Re from their OB v a lues if A b ecomes large. Surpr isingly, 
IWu and Libchabed (|l99lD and lZhang et all (|l997t) found 
that Nu is remarkably insensitive to NOB effects. To our 
knowledge NOB effects on Re have yet to be observed 
experimentally. 

In order to assess the va l idity of the OB approximation, 
iNiemela and Sreenivasa ^llooli) suggested that the fol- 
lowing three measures are basically equivalent: Busse's 
weighted sum of fractional deviations {Xt — Xh)/Xm for 
the relevant ma terial properties X, t he ra tio x = At/ At 
(introduced by IWu and Libchaben (|l99l[ )) of the tem- 
perature drops across the bottom and top BLs, and the 
relative change o f the density (3 A. The las t criter ion is 
the simplest, and INiemela and Sreenivasa 3 (l2003aD sug- 
gested on empirical grounds that OB conditions can rea- 
sonably be expected to prevail when pA is less than 
about 0.1 to 0.2. We will see below that the situation 
is more complicated and that there is a plethora of dif- 
ferent NOB effects. 

NOB effects for RB convectio n in water we r e meas ured 
and analyzed systematically bv lAhlers et al. 

Illooei- Of 

the relevant fluid properties, the kinematic viscosity i/ 
had the largest temperature dependence in this case. 
These authors compared Nu for three samples with dif- 
ferent L, but the same Ra, all with F = 1. The same Ra 
was realized in different ways: Either in a cell with larger 
L but A so small that NOB effects should be negligible, 
or in a cell with smaller L but larger A up to 40 K that 
had strong NOB effects. For the largest A the deviation 
of the center temperature Tc from the mean temperature 
Tm was onl y about 1.8 K (see fi g. [Tib), corresponding to 
X — 0.85. lAhlers et all ( 2006a[ ) calculated the increase 
of Tc theoretically by developing an extended Prandtl- 
Blasius BL theory, which takes the T-dependence of f 
and K into account. 

Surprisingly, for these water measurements again the 
NOB effects on Nu were tiny, see figure [TTk : less than 
1.4% even for A = 40K. To account for this finding, 
lAhlers et al] (j2006a[ ) derived an extension of the exact 
relation eq. ([5^ for the Nusselt number to the NOB case, 
proving the generally valid exact relation 



Nu 



2X1 



'oB KbAb + Kt At 



NuoB Xl[, + A^'., 



.A 



Fx -Fa. (34) 



1991 ) , and then using water and glycerol bv lZhang et al. 
19971 ) . These authors confirmed experimentally that the 



As can be seen, in each of the two factors Fx and Fa the 
symmetry-breaking different bottom and top BL proper- 
ties enter additively and thus tend to compensate each 
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to their temperature dependences. Thus the remaining 
NOB effects on Nu originate only from the quadratic and 
higher-order variations of the material properties with 
temperature. 

With the help of relation ^ lAhlers et al\ (|2006al ) 
could identify the origins of the NOB corrections for 
water: For this fluid it is mainly the temperature de- 
pendence of the thermal diffusivity k that is responsible 
for the NOB correction of Nu, while the NOB correc- 
tion of the center temperature Tc has its main origin in 
the temperature dependence of the kinematic viscosity 
J/. Instead in glycerol, which displays an extreme de- 
pendence of v on the temperature, the NO B corrections 
for Tr are much larger, both in experiment ( Zhang et al\ 
( 1997f)l and in the two-dim ensional numerical simulations 
of lSugivama et all (|2007t ). 



Theoretical calculations of T„ from the extended 



Prandtl-Blasius BL theory of lAhlers 6^0!! (|2006a^ for 
various water-like fluids with hypothetical temperature 
dependences of the material parameters were in good 
agreement with the correspo n ding numerica l resul ts of 
ISugivama efail (|2007l . l2008f) . lAhlers et al] (l2007l) ap- 
plied th e extended Prandtl-Blasius theory of Ahlers et all 
( 2006al ) to NOB effects in gases, finding good agreement 
for Tc with experiments using pressurized ethane gas, in 
which the material properties strongly depend on tem- 
perature. For this case a decrease of Tc relative to Tm 
and an increase of Nu in comparison to the OB case was 
found. Note that both effects arc in the opposite direc- 
tions as compared to the NOB effects in water or glycerol. 

In experiments by lAhlers et 'all (I2008D the measure- 
ments using ethane were extended to the region near 
the critical point, see fig. |TT] (b) and (d). On first 
sight the results for Tc seem surprising. NOB effects in 
liquid-like ethane (p > pc where pc is the critical den- 
sity) caused an increase of Tc relative to T™, whereas in 
gas- like ethane {p < pc) NOB effects made Tc smaller 
than Tm- The physical reason for this qualitative dif- 
ference was found in the opposite sign of the bottom- 
top asymmetry of the buoyancy strength due to the 
opposite temperature dependence of l3{T) for the two 
cases. In the liquid-like case the buoyancy, proportional 
to /3(r) = /3(T„0 + (T - r„0/3'(T,„), is larger at the 
bottom and smaller at the top, supporting the upris- 
ing warmer bottom plumes more than the down-coming 
colder top plumes. This brings predominantly hotter ma- 
terial into the bulk, leading to Tc > Tm- For gas-like 
ethane the buoyancy is larger at the cooler top, favor- 
ing down-going cold over uprising warm plumes. This in 
turn brings more cooler material into the bulk, resulting 
in Tc < Tm- It is the sign of the slope /3' of f3 at Tn that 
is the relevant quantity for this type of NOB correction. 

Of course, the extended Prandtl-Blasius BL theory 
fails when the thermal expansion coefficient (3, which is 
disregarded in that theory by construction as it treats the 
temperature as a passive scalar, shows an extreme tem- 
perature dependence such as close to a criticial point. 



ture dependence of /3, which leads to broken buoyancy 
symmetry, as an additiona l cause of NOB effe cts on Tc 
was verified numerically bv lAhlers et 

Even though the two cases p < Pc and p > Pc discussed 
above had Tc displaced away from r,„ in opposite direc- 
tions, they both yielded an enhancement of Nu above 
the OB value, i.e., in the same direction despite opposite 
deviations of Tc from Tm- This insensitiveness from the 
sign of Tc — Tm can be understood from Eq. ([M)) . Phys- 
ically, the BLs act as two thermal resistances in series, 
since Nusselt number deviations only depend on the sum 
of both, and it does not matter much which of them is 
reduced and which one is enhanced. 

The flow organization due to NOB effe cts in water and 
other liquids was studied numerically bv lSugivama et~all 
(1200811 . again with 2D RB simulations. It was confirmed 
that buoyancy, in particular the temperature depen- 
dence of the thermal expansion coefficient /3, is the main 
origin of the NOB effects on the center-roll Reynolds- 
number, which roughly behaves like Re/ReoB ~ 
(/3(Te)//3( Tm))^^^- Rey i iolds-n umber measurements us- 
ing water (jAhlers et all ( 2006a[ )) were still unable to re- 
solve NOB effects within their resolution of one or two 
percent even though A was as large as 38 K. 



VIII. GLOBAL WIND DYNAMICS 



A. Experiment 



LSC 



For cylindrical samples with F = 1 the 
circulation plane breaks the rotational invari- 
ance of the cell. This leads to interesting dy- 
namics which includes osc ill ations of the circula- 
pl ane (ICastaing et all (1X9891) 



tion 



img 

(1 19961). ICioni all (1 199?) 
iHeslot et all (Il987l). iNi cmcla et al 



Ciliberto et all 

'van Doom et all (12000^ 



(|2004 [) . lOiu and Tond (l2001aibl. 



2OOOD. iResagk oLl ( 20061) . ISreenivasan et all (I2002D 




Qiu et al. 



et al. 



iTakeshita ad (|l996f ). IXi and Xial (|2007l ). [Xi et all 
( 20061 )). These are caused by a torsional mode in which 



the orientation of the upper half of the LSC undergoe s 
azimuthal os c illatio ns (jFunfschilling and AhlersI ( 20041 ). 
iResagk et all (l2006h ) that are out of phase wi t h thos e 
of the lower ha lf (iFunfschilling and AhlersI (|2004l ). 
iFunfschilling et all (|2008l )). This is illustrated in the 
left part of Fig. [T^l which shows the LSC azimuthal 
orientations (characterized by the angle) 6b at the 
vertical position z ~ —Lj'^.^ the azimuthal orientation 
at z = 0, and the azimuthal orientation at z = -l-L/4 
(the origin of the z-axis is taken here at the cell center). 
These orientations were determined by measuring the 
azimuthal temperature variation of the side wall at 
the three vertical positions. Even casual inspection 
shows that Qt and Oh oscillate, out of phase with each 
other and perhaps with random amplitudes, about 
Qq- Quantitative analysis using correlation functions 
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The origin of this twisting mode is not known. How- 
ever, important insight is gained from the probabihty- 
distribution functions j) shown in the right part of 
Fig. [121 Here J is the probabihty of the azimuthal 
displacement 6*^ ^ = 9b.t — Oq, relative to 6*0, of 9b or dt. If 
the mode had its origin in a Hopf bifm'cation, then one 
would expect it to have a characteristic finite amplitude 
A which would lead to two peaks of p{0',,f) at or near 
9' = zLA. In contrast to this, the experimental result 
for p{0'fj f) is a near-Gaussian distribution with only one 
maximum centered at 0^ ^ = 0. Such a distribution is 
characteristic of a damped oscillator driven by a broa d- 
band noise source (see, for instance. iGittermanI ( 2005f )). 
This driving is attributed to the action of the small-scale 
turbulent fluctuations on the large-scale flow. 

The twisting mode does not influence the LSC ori- 
entation 6*0 at z = because there it has a node. 
However, the LSC breaks the rotational symmetry of a 
cylindrical sample and its circulation plane must some- 
how choose an azimuthal orientation. This orientation 
has bee n found to undergo sponta ne ous diffusive mean - 
dering (iBrown and Ahlerg (|2006alR l . ISun et all (|2005c[) . 
IXi et al. I 2006[ )) in the sense that its mean square az- 
imuthal displacement is proportional to the elapsed time. 
This diffusion is illustrated in the top of Fig. [13] by a 
long time series of 9o{t). It is also attributed to the ac- 
tion of the small-scale fluctuations on the large-scale flow. 
Measurements of its diffusivity have yielded results for 
the int ensity of the fluctuating force ( Brown and Ahler^ 
(l2006al lbh). 

Very recentlv IxTet al\ (|2008l ) and lZhou etai\ (|2008bl ) 
discovered a "sloshing" mode in a F = 1 cylindrical sam- 
ple that occurs in addition to the torsional mode and 
that was missed by the analy sis m ethod employed by 
iFiinfschilling and AhlersI (|2004[) and iFunfschilling et al\ 
( 20081 ) . This mode consists of an in-phase horizontal dis- 
placement of the entire LSC. The sloshing mode is said to 
be responsible for the oscillations seen in some of the lo- 
cal temperature measurements that had been attributed 
by some authors to periodic plume emission from the 
thermal boundary layers. 

The LSC also under goes re-orientation s both b y 
azimuthal rotations ( Brown and AhlersI (l2006d) . 
ICioni et al\ ( 1997t) ). and by cessations in which the 
LSC slows t o a stop and rest a rts in a random new 
orienta ti on (iBrown and AhlersI (|2006H) . 1 rown et al. 
(|2005bf ). IXi and Xial (|2007l )). A cessation is illustrated 
in the lower part of Fig. [TS] The left part reveals a 
sudden change in orientation. The right one shows that 
the temperature amplitude 6, and thus presumably the 
LSC itself, vanishes or comes close to zero at the time 
of the cessation. The angular change A6'c during a 
cessation has a nearly uniform probability distribution, 
indicating that, once the flow stops, it re-organizes itself 
in a random new orientation. The angular change 
of rotations has a power-law distribution, with small 
changes more likely then large ones. The time interval 



rotations is Poisson distributed, showing again that 
successive events are uncorrelatcd. 

On longer time scales, Earth's Coriolis force (at a 
latitude of 34° in the Northern hemisphere) was found 
to cause a net clockwise (as seen from above) rotation 
of the LSC orientation on average once every 3 days, 
and to align the LSC in a preferred orientation close 
to West. The measured probabihty distribution ^(6*0) 
agreed quantitatively with a model calculation of the 
Coriolis-f orce interaction tha t involv ed no adjustable pa- 
rameters (Brow^^nT^hleri (20063) )■ The net clockwise 
rotation also was consistent with this model. 

The LSC in other geometries and aspect ratios adds to 
the richness of this phenomenon but is beyond the scope 
of this review. We mention briefly that DNS revealed two 
counter-rotating n ear-circular rolls, stacked vert ically one 
above the other ( Verzicco and CamussH (|2003[ )) as the 
dominating feature of the LSC in cylindrical samples 
with F < 0. 5 . Rece nt experimental evidence provided by 
IXi and Xial (|2008bl ) for F = 0.5 and 0.33 indicates ran- 
dom temporal successions of one-roll and two-roll states, 
with the two -roll state becomin g more prevalent as F 
is decreased. IXi and Xial (|2008aD found a strong aspect 
ratio dependence when comparing the azimutal motion, 
reorientation, cessation, and reversal of the large-scale 
circulation in cylindrical RB samples with F = 2.3, 1, 
and 0.5. In samples of square or rectangular cross sec- 
tion the LS C is locked in a p r edom inant orientation along 
a diagonal ( Dava and Eckel ( 200ll )). This is now under- 
stood in terms of the pressure gradients that arise when 
the rotational invar i ance of a c ylindrical sa n iple is broken 
([Brown and AhlersI (|2008a[) ). Izhou et al\ (|2007a[) stud- 



ied temperature and velocity oscillations in a rectangu- 
lar cell, too, finding that the temperature oscillations 
scale differently with Ra than the velocity oscillations, 
which are affected by the cell geometry. Of great in- 
terest, but largely unexplored, is the fiow structure that 
will be found in large-F systems where numerous con- 
vection rolls potentially can co-exist next to each other. 
Recent papers addres s ing th is issue numerically include 
that bv lHartlep et al\ I 2005t l. Recent experimental work 
going in this direction i i icludes the measurements by 
iNiemela and Sreenivasa 3 (l2006aD . 



] B. Models 

Stocha stic models of fiow rev ersal ha ve been pro - 
posed bv lSreenivasan et~ai\ ( 2002h . and bv lBenzil ( 20051 ). 
They treated diffusion of the LSC strength in a poten- 
tial well, but there was no physical motivation for the 
shape of the potential that was used and the model pa- 
rameters were chosen phenomenologically. They also 
did not address the azimuthal dynamics of the LSC. 
Two other models describe the LSC with determin- 
istic ordi nary differential equatio n s tha t have chaotic 
solutio ns ( Fontenele Araujo et al\ ( 20051 ). iResagk et al\ 
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guably relevant aspects of the Navier-Stokes equations 
and making various approxi mations. The model of 
iFontenele Arauio et al\ ( 2005[ ) is based on assumptions 



about the life times Tp of plumes. It is estimated that 
plumes with a sufficiently large Tp will be carried over by 
the LSC to the far side where their buoyancy tends to 
act against the prevailing flow. This physical process can 
lead to cessations. Since this model does not contain an 
azimuthal mode, it can not describe the rich azimuthal 
dynamics of the physical system. 

Recently a model consisting of two stochastic ordi- 
nary differential equations, one each for the circulation 
strength U and the azimut hal circulation-plan e orien- 
tation 9o, was developed by iBrown and Ahlerd (l2007al 
l2008bt ). Starting with the Navier-Stokes (NS) equation, 
it is argued that U is driven by the buoyancy term and 
hindered by the dissipation in the viscous boundary lay- 
ers near the walls. A phenomenological stochastic driv- 
ing term representing the interaction between the small- 
scale turbulent fluctuations and the large-scale circula- 
tion is added to the deterministic model. For Oq the only 
driving is the turbulent fluctuations, but the component 
of the nonlinear term in the NS equation that describes 
the rotational inertia of the LSC provides damping and 
couples the two equations. Assuming that the experi- 
mentally accessible azimuthal temperature amplitude 5 
is proportional to U near the side wall, the model be- 
comes 
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The intensities of the stochastic forces fs{t) and fg{t) 
are obtained from model-independent measurements of 
the diffusivities of S and ^o- 

The potential corresponding to the 6 equation has a 
(5^/^ rather than the usual quadratic nonlinearity. How- 
ever, this does not change its qualitative structure, which 
has an unstable fixed point at S = and a stable one at 
6 ~ Sq. The driving will cause diffusion in the vicin- 
ity of the stable fixed point, with occasional excursions 
to 6 = corresponding to cessations. During a cessa- 
tion (when 6 is small) the damping term on the rhs of 
the 00 equation, and thus the angular momentum of the 
LSC, become small, and it becomes easy for the stochas- 
tic driving to cause relatively large angular changes AOc- 
The model agrees well with many experimental r esults 
for the LSC dynamics (jBrown and AhlersI (|2008bh ). in- 
cluding a time interval between cessations of 1 to 2 days, 
a near-uniform distribution p{A6c) for cessations, and 
the dependence of Sq on Re and Ra expressed by the re- 
lation (So /A) X (Ra/Pr) cx Re^^^. There is, however, a 
significant disagreement between the model and the mea- 



for 6 < Sq/2. It has been suggested ( Brown and Ahle"ri 
(j2008bl )) that this is caused by the neglect of thermal 
conduction across the thermal BLs which is expected to 
become important when S is small. 

Very recently the abo ve model was extended by 
IBrown and Ahle^ 12008^ by including various pertur- 
bations that break the rotational invariance of the sam- 
ple. These include (a) the effect of Earth's Coriolis force, 
(b) an elliptic rather than circular horizontal cross sec- 
tion of the sample, (c) a tilt of the sample axis rela- 
tive to gravity, and (d) a small horizontal temperature 
gradient at the top or bottom plate. It turns out that 
(a) and (b) only influence the 9q equation and not the 
S equation. Perturbations of this type, although they 
introduce a preferred orientation of the LSC, leave the 
Reynolds number and the frequency of cessations largely 
unchanged. However, perturbations like (c) and (d) 
which affect 6 influence Re and suppress cessations. A 
tilt or elliptic eccentricity of sufficient magnitude creates 
a new oscillatory mode of the LSC which is different from 
the torsional oscillation in that the phase of the oscilla- 
tion is uniform along the height of the sample. For the 
tilted sample this mode wa s found in recent experiments 
([Brown and AhlersI (|2008at )). 



IX. ISSUES FOR FUTURE RESEARCH 

As seen in this review, since Siggia's 1994 article in An- 
nual Reviews of Fluid Mechanics tremendous progress of 
the understanding of the turbulent Rayleigh-Bcnard sys- 
tem has been achieved by experiment, theory, and nu- 
merical simulation. However, it has also become clear 
that our understanding is far from complete. In the fol- 
lowing we summarize what we consider as major issues 
for future research. 

Presumably the most important challenge is to clarify 
whether and, if so, where in parameter space the ultimate 
state of convection exists. Estimates suggest that (for Pr 
near 1) a transition to such an ultimate state should oc- 
cur at Rayleigh numbers around 10^'^ to 10^^. Beyond 
that transition the Nusselt number should increase more 
rapidly with Ra than below it. Though the Grenoble 
experiments suggest such a transition near Ra = 10^^, 
neither the Oregon/Trieste experiments nor numerical 
simulations do so. The reason for the discrepancy is un- 
resolved at present. It is of utmost importance to clarify 
this issue in order to allow extrapolations of the heat 
transfer to the very large Ra regime of geo-physical and 
astro-physical interest. Perhaps related to this issue is 
whether a coherent LSC continues to exist at very large 
Ra, or whether it is totally overwhelmed by fluctuations. 
Several experimental efforts are planned to try to answer 
these question. 

Rather than focusing on global quantities such as the 
heat flux, which may be very difficult to control in the 
large-scale setups necessary in an ultimate regime, a 
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detailed analysis of the top and bottom BLs, whose struc- 
ture should reflect such a transition. Here the challenge 
is the opposite: Extremely small structures must be spa- 
tially resolved. Even for an about seven-meter high sam- 
ple (such as the Barrel of Ilmenau) the thermal boundary 
layer is only about 1mm thick when Ra = 10^^. Helpful 
strategies may be to try to trigger the transition in the 
BL, by controlled roughness or even by moving parts. 

Another major question to be clarified is the role of the 
plate corrections: Though major efforts have been under- 
taken to correct for the finite plate heat-conductivities, 
beyond Ra ~ 10^ there still seems to be a discrepancy 
of 20% or more between experimental measurements of 
Nu and numerical data for constant plate-temperature 
boundary-conditions, see figure [SI In fact, this question 
may even be related to above one on the existence of the 
ultimate regime. 

Also the three-dimensional dynamics of the large scale 
convection roll will need further attention and analysis. 
As we have seen in section lVHIl it is rather rich, including 
torsional modes, rotation, cessation, and sloshing. Here 
the key question is: What features of the LSC can be 
described through a deterministic model somehow based 
on the Oberbeck-Boussinesq dynamics or some force bal- 
ance and what features need stochastic elements for a 
description? Routes for further research on this subject 
will include experiments and numerical simulations with 
modified geometries as e.g. cylindrical samples with a 
horizontal axis. 

Another issue for further research is the exploration 
of the aspect-ratio dependence of the flow. Though the 
F-dependence of Nu is found experimentally to be re- 
markably weak, a full understanding of the RB system 
requires further theoretical efforts also in that direction. 
For the characterization of the corresponding flow organi- 
zations detailed PIV measurements could play a crucial 
role. These will also shed light on the Ra and Pr de- 
pendences of the various Reynolds numbers that one can 
define, and their connections. 

From the theoretical viewpoint, the GL theory has 
provided a good guideline for the understanding of 
Nu(Ra,Pr) and Re(Ra,Pr) and even allowed various pre- 
dictions, but clearly the theory has its limitations: 

• It builds on the Prandtl-Blasius BL theory for the 
temperature as a passive scalar. The buoyancy 
term is skipped by construction in the Prandtl BL 
equations, but obviously it is of high importance for 
thermally driven fiow, leading to plume detachment 
from the thermal BLs. This detachment mecha- 
nism as a time dependent BL separation process 
needs further study and analysis of its parameter 
dependences. The GL theory only represents the 
global effect of the plumes and the resulting self- 
organizing buoyancy-driven flow, namely the large 
scale wind. 

• The three-dimensional (torsional) dynamics of the 



the GL theory (or any other model). Although 
it was established experimentally that this oscil- 
lation is a stochastically driven mode, its origin is 
presently unknown. 

• Experiments suggest that Nu and Re are not as 
intimately related as the relations (|2ip and of 
the GL theory suggest. For instance, by tilting the 
system Re can be enhanced considerably whereas 
Nu changes very little. 

• Figure 3] implies that there are more sudden tran- 
sitions between the various regimes in reality than 
within the GL theory. 

• Finally, in the Nu oc Ra^^'^ regime the experimen- 
tally found co-existence of a measured Reynolds 
number oc Ra^^^, achieved in Santa Barbara, and 
the GL prediction Re cx Ra^^^ for the global wind 
has to be clarified and resolved. 

More theoretical work is necessary to improve on the GL 
theory. This unfortunately may only be possible by sac- 
rificing its conceptual appeal and simplicity. 

An important further extension of the present theo- 
retical understanding of RB convection, but of course 
also of the numerical simulations and the experiments, 
is the extension towards rotating RB convection. With 
rotation, the (Ra,Pr)-parameter-space is extended to 
(Ra,Pr,Ro). where the Rossby number Ro is defined as 
Ro = {2Q)^^ ^/ PgA/ L, i.e., as ratio of the time scale 
given by the rotation rate CI and the one given by buoy- 
ancy. The obvious questions to address are: What is the 
dependence of the Nusselt number on the control param- 
eters Ra, Pr, and Ro, i.e., what is Nu(Ra,Pr,Ro)? And 
how do the large scale convection roll and the Reynolds 
number react to the rotation, i.e., what is Re(Ra,Pr;Ro)? 
When will the large scale convection roll break down? 
How are the top, bottom, and side-wall kinetic and ther- 
mal BLs modified through the rotation? Only a small 
fraction of the parameter space (Ra,Pr,Ro) has hith- 
erto been explored. Given that rotating turbulence is 
known for its counter-intuitive features (e.g. the Taylor- 
Proudman effect), we expect many surprising results 
ahead of us. We do not want to give an extensive lit- 
erature review on rot ating RB here, bu t good starting 
points ar e the book by Greenspaiil ( 1990t ) or the classical 
article bv iRossbvl ( 19691 ) himself. 

Future developments undoubtedly will also include the 
extension of the current state of the art of RB convection 
to complex fiuids, to fluids containing bubbles or suspen- 
sions of particles, to fluids undergoing phase transitions, 
and so on. But clearly these topics are beyond the scope 
of the present review. 

Fifteen years after Siggia's Annual Review of Fluid Dy- 
namics article on RB convection, and inspite of the huge 
progress achieved during this time, we still close with a 
similar statement as he did: The current rate of experi- 
mental and numerical advances will again soon antiquate 
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FIG. 3 U pper: Phase diaRram in the Ra — Pr plane ac - 
cording to iGrossmann and Lohsd (|200d . 1200 ll . |2002| . 12004 ): 
The upper soUd Une means Re = Rbc, the lower nearly par- 
allel solid line corresponds to tu,BL = tu,buik\ the curved 
solid line is eg^BL ~ ee,buik', and along the long-dashed line 
All = Xe, i.e., 2aNu = V Re. The dotted line indicates 
where the laminar kinetic BL is expected to become turbu- 
lent, based on a critical shear Reynolds num b er J^e * ~ 420 
of the kinetic BL, cf. iLandau and Lifshitj (| 19871 ). Data 
points where Nu has been measured or numerically calcu- 
lated have been included (for asp e ct rat ios F ~ 1/2 to 1): 
Chavanne et al\ (Il997f): cy a n diam onds : 
red circles: iNiemela et al\ (l2000al): blue 



green squa res: 



ICioni et al\ (Il997 



stars: lAhlc rs and Xul (|200ll); re d stars: iFunfschilling et al\ 
(120051). [Nikolacnko et al\ (|2005l ): bro wn triangle s down : 
IXia et al\ (|2002l ); cyan t riangles down: ISun et al\ (|2005bl '): 
maroon tria ngles right: Idu Puits et al\ ( 2007bl ) :~~violet tri- 
angles up: IVerzicco and Camussi~(ll99!: ) (numerica l sim- 
ulations) ; blue squares: iKerr and Herrin"3 (l2000t) (nu- 
meric al simulations) ; magenta t r iangle s up: lAmati et al\ 
l|2005l ). IVerzicco and SreenivasanI (|2008l ) (numerical simula- 
tions). Note that some of the large Ra data probably are 
influenced by NOB effects. Lower: an enlargement of the 
upper figure. 
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FIG. 4 Reduced Nu for F = 1, obtained using water (Pr ~ 
4.4) and copper plates, as a function of Ra. Open symbols: 
uncorrected data. Solid symbols: after cor rection fo r the fi - 
nite plat e conductivity. C ircles: iFunfsclii lline: et al\ (|2005l ). 
Squares: ISun et al\ (j2005al ). The downwards and the upwards 
triangles are upper and lower bounds on the actual Nusselt 
number at large Ra; the diamonds originate from an estimate, 
see the text for details. The solid line is the GL prediction 
(IGrossmann and Lohsd (|200ll) ). 
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FIG. 5 Top figure: Nu( Pr) for Ra = 6 ■ 10^ and T = 1 from 
numerical simulations bv lVerzicco and Camussil (Il99£ ) (black 
circles), from experiments with mercury by " iRossbvl (|l969t ) 
(blue diamond), and from the experiments with sodium by 
iHoranvi et al\ (| 19991 ) (red square) . Th e straight solid line is a 
fit to the numerical data with Pr < 1 (|Verzicco and Camussil 
(| 19991 )). giving Nu = 8.1Pr°-^*±° °^ The exponent is in 
agreement with the low-Pr expectation 1/8 of the GL the- 
ory. The gr een upper triangle s are t he numerical data for 
Ra = 10^ by iKerr and Herrind (|2000l ) the dark- green circle 
results from the experimental data of ICioni et al\ (|l997l ) for 
the same Ra = 10^, a nd the violet d iamon ds are numerical 
results for Ra = 10*^ bv lBreuer et al\ (|2004l ). The three thick 
red lines are the results from the GL theory eqs. (|21|l and (|22p 
for the three Rayleigh numbers of the numerical data sets, 
namely Ra = 6- 10^, Ra = 10®, and Ra = 10^, bottom to top. 
Bottom figure: The reduced Nusselt number NuRa^^''* as a 
function of the Prandtl number for the two Rayleigh num- 
bers 1.78 X 10^ (upper set) an d 1.78 x 10^ ( lower set) in the 
large-Pr regime. Open ci rcles: lAhlers and Xu ( 200l| ). Solid 
symb ols: IXia et al\ (|2002| ) . Various organic fluids were used. 
From IXia et al\ l|2002l ). 
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FIG. 6 Top figure: Nu/Ra^«/^ with 7^/7 = 0. 323 as a func- 
tion of Ra. Black soUd circles: F — 0.5 (|Niemela et al\ 
2000al')^ after a correction for p late and side-wall effects 



..Niemela and SreenivasanI (l2006bl ll. Black open circles: F = 
0.5 ( IChavanne et al\ 1120011)'). Blue solid squares: F = 1 



([Niemela and SreenivasanI (l2003a|lV Purpl e solid diamonds: 
F = 0.67 and 0.43 (iNikolaenko eta!] (120051')'). Open s quares : 
1 ^ r < 3, Pr = 0.7 (|Fleischer and GoldsteinI (120021)'). Red 
stars: numerical results for F — 1/2, Pr = 0.7 ( Amati et al\ 
(|2005l )) for constant-temperature boundary conditions at the 
plates. Blue stars: numer i cal re sults for F = 1/2, Pr — 0.7 
(IVerzicco and SreenivasanI (|2008l )) for constant-heat-flux BCs 
at the lower plate and constant-temperature BCs at the upper 
plate. Dotted red line: a power law with ysff = 1/3 fitted to 
the red stars. Solid (dotted) black line: The GL prediction 
for Pr = 0.8 (Pr = 29). Bottom figure: The Prandtl numbers 
corresponding to the data in the top the figure. In addition 
Pr for the measurements for F = 1 with water in Fig. [J] are 
shown as black diamonds. The dashed and dotted lines in the 
bottom figure are estimates of the location of the transition 
to the Kraichnan regime for F = 1, assuming critical bound- 
ary layer Reynolds-numbers Re* — 220 and 440 respectively 
(since the parameters of the GL theory have been determined 
only for F = 1, a prediction of Ra* for smaller F unfortunately 
is not available). 
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FIG. 7 (a): Normalized horizontal velocity u/ umax as a func- 
tion of the vertical position pz = z/L (from lOiu and Tond 
(|2001bf )'). Note that the data extend only up to ~ 0.85; 
thus the minimum near the sample top is not shown, (b): Re- 
duced Reynolds numbers Re"/Ra^'''^ (circles) and Ke'^ /Ka^^^ 
(solid triangles), defi ned in the text, as functions of Ra (from 
iQiu and TonS (|2002l )). 
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FIG. 8 (a) and (b): 
(a): Re"' for Pr = 
triansles) , and Pr 



Brown et al 



Brown et al 



Reduced Reynolds numbers Ke/Ka}^^ . 
4.38 (circles), Pr — 5.55 (up-pointing 
= 3.32 (down-pointing triangles) from 
Stars: Re" for Pr = 28.9 from 
Plusses; Re" from lQiu and Tond (j2002h . 
Pr ~ 5.4. Dashed lines (from top to bottom): GL predictions 
for Pr = 3.32, 4.38, 5.55, and 28.9. (b): Re"' (solid symbols) 
and Re" (open symbols) for Pr — 4.38 (circles), Pr = 5.68 
(up-pointing triangles) and Pr = 3.26 (down-pointing trian- 
gles). Dashed lines (from top to bottom): GL predictions for 
Pr = 3.26, 4.38, 5.68. (c): The ratio of Re to the GL pre- 
diction as a function of Pr for Ra < 3 x 10®. Open circles: 
Re"', Pr = 4.38. Stars: Re", Pr = 29.8. Squares: Re" from 
iLam et al\ (|2002l '). 5.6 < Pr < 1206. 




FIG. 9 (a) Local, only time-averaged, temperature and 
temperature-fluctuation proflles at Ra = 7.7 • 10^^ and F = 
1.13 as functions of the height z (normalized by Xf), measured 
below the center of the top plate. The fluid is air at atmo- 
spheric pressure, the cell height is 6.3 m. From ldu Puits et al\ 
(|2007bD . (b) Rescaled mean velocity profiles U{z) above the 
center of the plate for Ra = 1.25 • 10^ (shortest data set) up 
to Ra = 2.02 • 10^° (longest data set) for RB fiow in water. 
u is normalized by Umax and z by the kin etic BL thickness 
Xu{*). Figure taken from lSun et al\ (|2008l '). 
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FIG. 10 Scaling of the kinetic BL thicknesses X^i*) and 
A™(*) as functions of (a) Ra and of (b) Re, taken ab ove the 
center of the lower plate. Taken from lSun et al\ (|2008l ). In (a) 
also the scaling of the thermal BL thickness A|'(*) is shown. 
It has been deduced from the temperature profile above the 
center of the lower plate. The lines are power-law fits. 
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FIG. 11 (a) and (b): The ratio Nu/Nuqb as a function of 
A for fixed Tm- (c) and (d): Tc — Tm as a function of A, 
also fixed Tm- (a) and (c): Water at Tm = 40.0° C, T = 1.0. 
Solid circles: L = 49.7 cm. Open circles: L — 24.8 cm. Solid 
squares: L = 9.2 cm. (b) and (d): ethane at a pressure of 
55.17 bar for F = 0.5. Open symbols: = 35.0° C below the 
critical isochore temperature = 38.06° C at this pressure. 
Solid symbols: T™ = 41.0°C, above T^. 




FIG. 12 Left: The LSC orientation as a function of time at 
z = — L/4 {9b, dotted line), z = {9o, triangles), and z = L/4 
{9t, solid line). Right: The probability distribution p{9' /2n) 
as a function of 6' /2tt. Here 6' — 9t — 9o (open circles) or 
9' = 9b — 9o (s olid circles) . Solid line: a G aussian distribution. 
Adapted from lFunfschilling et al. \ (|2008D . 
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FIG. 13 Top: A time series of 9o ( t) span ninK about 11 days 
(adapted from lBrown and" AhlersI (|2006bD '). Bottom left: the 
LSC orientation 9o{t), bottom right: the LSC azimuthal tem- 
perature amphtude 5(t) at the side wall, both during a cessa- 
tion (adapted from iBrown et al\ (|2005bf )'). 



